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ABSTRACT 


A method has Been developed for prediction of azeotropic 
data at different pressures from knowledge of the azeotropic 
data at a particular pressure. The C-auss-Seidel iterative 
technique is used on computer for solving the mathematical 
equations developed. The non ideality in the vapor phase is 
incorporated m terms of the vinal equation of sbatc and the 
excess free energy of mixing is calculated using any one of 
the correlations such ass R-K equation, Van-Laar equation and 
Scatchard-Hammer equation. The method has Been verified using 
data from the literature for the systems s Ethyl alcoho 1-Ethyl 
acetate, Benzene-Cyclohexane, Hexane-Benzene, Methyl elcohol- 
Benzene, Ethyl alcohol-Methylcyclohexane and Ethyl alcohol-water. 
The predicted azeotropic temperatures and compositions are found 
to he m very good agreement with the experimental values except 
for highly nonideal systems comprising polar components in which 
maximum absolute error is 4 - 1 °, 



CHAPTER 1 


INTRODUCTION 

An azeotrope is a constant toiling mixture in which 
the vapour and liquid phases in equilibrium have the same 
composition, because of which components can not be separated 
by simple distillation. 

Pormation of azeotropes is a common phenomenon in 
industrial separation processes which makes one interested 
in the study of their vapour-liquid equilibrium behaviour. 

It is often observed for certain systems to exhibit either 
vapour pressure extremum under isothermal conditions or 
boiling temperature extremum under isobaric conditions. The 
azeotropic composition corresponds to either of these extre- 
mum points. 

Azeotropes may be classified broadly into Homo and 
He tero azeotropes depending upon the number of liquid phases 
present m the azeotrope. In the case of Homo azeotropes, 
components of the mixture are completely miscible with each 
other resulting m a single liouid phase whereas m the case 
of Hetero azeotropes, components are partially miscible resulting 
in more than one liquid phase. The mixtures however may change 
from one form to the other on varying the temperature over a 
certain range. Azeotropes may be further classified into 
positive and negative types* A positive azeotrope has a minimum 
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boiling point at constant pressure or maximum vapour pressure 
at constant temperature. Conversely, a negative azeotrope 
lias a maximum boiling point and a minimum vapour pressure as 
shown m Pig. (l-l). Saddle (positive-negative) azeotropes 
are those azeotropes m which both maximum and minimum boiling 
points appear for a particular system at different compositions. 
They usually appear m multicomponent systems. 

The azeotrope formation generally occurs when boiling 
points of the components lie within a certain range. Azeotropes 
are formed because of high nonideality m mixtures. If compo- 
nents differ largely in their molecular structures, sizes and 
chemical nature, the -chances of formation of azeotropes are 
quite significant. It is also observed that m the case of binary 
mixtures, if the boiling points" of the two components are iden- 
tical, the azeotropic composition is approximately 50 mole per 
cent and shifts to the right or left when the boiling point of the 
the component 2 is below or above that of 1. Azeotrope forma- 
tion over a pressure range will depend upon the temperature 
coefficients of the vapour pressures, i.e. dP/dt of the two 
components. If the coefficient for component 1 (more volatile) 
is much greater or less than that of 2, the azeotrope dis- 
appears before the critical region is reached as shown m 
Pigure (1-2). The dotted lines show the locus of the azeotropes. 
If the coefficients of the components are about the same, the 
azeotrope is terminated m the critical region as shown in 
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Figure (1-3). This is called the critical azeotrope. In 
the critical region the critical locus curve (i.e. T' c vs x) 
is characterized by a minimum temnerature point which not 
necessarily represents the critical azeorrope. A point of 
maximum and or minimum pressure in the critical locus curve 
(Pc vs x) may also exist when the composition of the critical 
azeotrope is rich m one of the components. Azeotropes formed 
m the critical region, only for those compounds whose cri- 
tical temperatures do not differ by more than a characteristic 
value depending upon the chemical nature of the compounds under 
study. 

The azeotropic behaviour of a system is studied by 
knowing its pressure-temperature and composition relationship 
or m other words the azeotropic P-I-x space. Once the 
azeotropic P-T-x space is known, one can separate the components 
either by varying the pressure or temperature, or by changing 
the composition. The composition can be changed by adding 
one of the components m excess or a third component to the 
mixture. 

Most of the azeotropic pressure temperature composition 
data available m the literature are at atmospheric pressure. 
Very limited data are available at pressures other than 
atmospheric pressure. In the absence of a good method or a 
correlation for predicting the azeotropic P-T-x space over 
a wide range of pressures an attempt has been made m this 




study to develop a theoretically sound method for accurate 
prediction of these data. 



CHAPTER 2 


LITERATURE REVIEW 


Dalton was the first person to observe the occurrence 
of a minimum vapour pressure while distilling aqueous solutions 
of hydrochloric acid and nitric acid. He observed that the 
boiling temperature and composition of distillates remained 
constant and also the boiling temperature was always higher 
than that of the higher boiling component. Roscoe attributed 
this phenomenon to the formation of a new chemical compound 
which caused the lowering of total vapour pressures. Between 
1881 and 1884 the Russian chemist Konovalov observed the occurre 
nee of minimum and maximum vapour pressures in the mixtures for 
which equilibrium liquid and vapour compositions we ro identical, 
and the maximum vapour pressure associated with the minimum 
boiling temperature and vice-versa. Later Storonkm observed 
this phenomenon m multicomponent mixtures. In 1911 Wade and 
Merriman introduced the term AZEOTROPE to designate all such 
mixtures irrespective of the number of components. Azeotrope is 
a Greek word which means 'difficult to separate'. 

Since Azeotropism is one of the most common manifestation 
of nonideality m liquid solutions, the pioneer thermodynami- 
cists like Duhem, Marqules, Lehfeld and others attempted to 
explain the occurrence of azeotropes theoretically. 
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In 1910 the Russian scientist VREVSICI formulated a 
general rule for the effect of pressure on azeotropic compositic 
of a binary system. According to this rule, for a positive 
azeotrope the concentration of the component with higher 
evaporation enthalpy is increased with pressure, on the other 
hand, for negative azeotropes under isothermal conditions the 
change m pressure produces a decrease in the concentration 
of the component having higher evaporation enthalpy. This 
rule can be thermodynamically explained by equation (2.1) whic v 
has boon derived with two assumptions, namely (i) Molar volume 
of the liquid is assumed to be negligible m comparison with 
the molal volume of the vapour, and ( 11 ) Vapour phase is 
assumed to behave like an ideal gas. 


dx.j 

dlnP = 








( 2 . 1 ) 


where? .All = x^ Z + ^2 ^ ^2 an ^s ~~ ^ and Z ^ are 

the partial molal vapourization enthalpy of the components m 
the mixture. Ti/henever activity coefficient is greater than 
unity, the denominator of equation (2.1) is found to be 
negative and hence according to the equation (2.1) The mole 
fraction of component- 1 m the azeotropes increases which 
pressure provided /IH-. is greater than 

Lecat (1918-1928) carried out numerous experiments on 



azeotropes using highly pure chemicals and compiled extensive 
data known till then. Later in 1947 and 1948 Horsley compiled 
all experimental data known to date on "binary and ternary 
zeotropes (non-azeotropes) as well as azeotropes. 

In 1929 Lecat showed that the relation "between the 
pressure and "boiling temperatures of the azeotrope is similar 
to that for a pure compound. It may be represented to a close 
approximation by the equation (2.2), 

log P = A - B/T (2.2) 


where A and B are constants characteristic of the system. 

21 

Prigogme employed claussius-Clapeyron equation (2.3) +- 
predict pressure-temperature relationship of azeotropes assum 
heat of vaporization of the azeotrope to be independent of 
temperature i.e. 


dT _ 
dl “ 



(2.3) 


where, = heat of vaporization of tho azeotrope and = 

volume change m vaporization of the azeotrope. 

Later Nutting and Horsley' made use of the Cox-chart 
to predict azeotropic pressure-temperature relationship. 

An approximate method for estimating the effect of 
temperature on azeotropic composition was suggested by Carlson 
and Colburn . This method was based on the assumption that tho 
ratio of the activity coefficients is independent of temperature 
and also that the vapour phase to behave like an ideal gas. 
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Hence the ratio of the activity coefficients is given by 
equation (2.4), 



where, , P^,;, ,,y are the activity coefficients and vapour 

pressures of components 1 and 2. According to this method 
y-i 

values of -rp- are plotted against the liquid composition, x, 

) 2 

for a known temperature and on the same plot, values of I^/P^ 
are plotted against temperature using the same ordinate scale 
(Figure 2.1). From this plot azeotropic composition at any 
de sired temperature can be obtained by making use of the 
azeotropic condition i.e. ^/} ) 2 = P 2 /P^» Tlie ^thod 1S 
demonstrated by the dotted lines in Figure 2.1. For estimate 
of the total pressure corresponding to various azeotropic 
compositions and temperatures the following empirical expression 
was suggested i 


P’ (x 1 P 1 + x 2 P 2 ) 
z = ""(xp^ +"2j rp 


(2.5) 


where, P z , P^ = total vapour pressures of the azeotrope at the 
desired temperature, t and the reference temperature, t*. 

P.j, P£ = vapour pressures of the pure components 1 and 
2 at the reference temperature, t’. 

, x 2 = mole fractions of components 1 and 2 m the 
azeotrope at t. 


x.j,x^ = mole fractions of components 1 and 2 at t 1 . 




( Composition ) 


FIG- 2.1 COULSON COLBURN METHOD FOR PREDICTION 
OF AZEOTROPIC COMPOSITION . 
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However the equation (2.5) finds little use as it 
often renders high inaccuracy m estimation of the total 
pressure, 

Skolnic^ (1951) developed a graphical method for pre- 
dicting pressure-temperature relationship of the azeotropes 
using the nomographs of Lippincott and Lyman . These nomographs 
give vapour pressure-temperature relationship for pure compounds 
having normal boiling points between -50 and 550 °C. Prior to 
it a modified form of the integrated Claussius-Clapeyron 
equation ( 2 . 6 ) was developed to depict the vapour pressure- 
temperature relationship. 


A+ - (275.1 + t) (2,8808 log. 
4 VC _ 0 4 0.15 (2.6808 - log 




( 2 . 6 ) 


where, At = temperature m °C to be added to the observed 
tempera fcure. 

t = observed temperature 

(2.8808-log p) = Log of observed pressure substracted from 
Log of standard pressure (760 mm of Hg. ) 

0 = quantity proportional to entropy of vaporization 

at 760 mm of Hg. 

The solution of equation (2.6) is found to be tedius 


and time consuming. To overcome this, Lippincott classified 
all the compounds into several groups according to their 
resemblence m physical and structural properties and presented 
nomographs to get vapour pressure-temperature relationship 
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knowing only the normal boiling temperature of the compound 
under consideration. It was observed that the average absolute 
deviation of the boiling temperature of pure compounds obtained 
from the nomographs from the experimental values is 0.9°U. 

Use of Lippmcott ’ s nomographs was extended to the 
production of azeotropic pressure-temperature relationship. 
Because an azeotrope is a mixture of either two or more compo- 
nents, usually of widely different physical and structural 
properties, it is classified under component groups. Let us 
illustrate the method for an azeotrope whose components belong 
to groups 3 and 8. If the normal boiling point of the azeotro o 
is denoted by 9 P C and the reference pressure by P , Then the 
pseudo boiling points of the components at this pressure are 
given by t^ and ^ as shown in Figure 2.2, Then correspond inf 
to these temperatures points C and A are plotted in Figure 2,3 
and D corresponds to normal azeotropic point. Finally distance 
CA is divided at point B m the ratio of the compositions of 
the two components known at standard, pressure. Hence path 
DB represents pressure -temperature relationship for the azeo- 
tropic system. 

This is obviously an approximate method as it is known 
that composition changes significantly with pressure. However, 
it was observed that for most of the azeotropic systems, this 
approximation does not lead to significant error in estimation. 
The deviation of the predicted temperature from the experime,.' 
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value is of the order of 1°C, The method is useful only at 
subatmo spheric pressures and has limitations namely, the normal 
■boiling temperature of the azeotrope must he known and if the 
components of the azeotrope belongs to same group then this 
method can not he used for prediction of the P-T-x data. The 
method obviously fails for systems for which azeotropes form 
only at pressures above atmospheric. Once the pressure- 
temperature relationship is known, the azeotropic composition 
is obtained from the linear relationship assumed by Skolnic 
for low pressure s; 

log x = A - BT (2.7) 

where, x = mole percent of one component m the azeotrope 
T = boiling temperature in °C of the azeotrope 
A,B = constants for a particular system. 

This technique was also tried for some multicomponent 
system but yielded poor agreement compared to binary systems. 

Bandara^an^ of IIT/Madras also suggested a method for 
prediction of pressure- temperature relationship m Binary 
systems comprising non-polar as well as polar compounds from 
the pure component properties, such as molal latent hects of 
vaporization, molar volumes, vapour pressures etc. Heed’s 
method was used for prediction of the azeotropic behaviour of 
non-pSlar mixtures. However, when the method was used for polcir 
mixtures, the deviation of the predicted values from the 
experimental ones were quite large. Hence a method was developed 
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to relate these deviations to the dipole moments of the 
components. The agreement m boiling temperatures was found 
to range from within + 0.3 to 3°C. 

Thus an extensive literature survey on the existing 
methods for prediction of azeotropic pressure-temperature- 
composition data reveals that the methods are either empirical 
m nature or very complicated and time consuming. Most of the 
methods have been used at the low pressure range and yield 
inaccurate results. Hence it is decided to develop a method, 
for predicting the azeotropic P-T-x from the knowledge of one 
azeotropic point, which is theoretically sound having least 
assumptions and yet would yield better agreement compared to 
existing methods. 


#•*** 



CHAPTER 3 


THEORETICAL CON SI DE -ATMS 

According to the phase rule there is only one degree 
of freedom. for a "binary azeotrope. Eence the P-T-x data for 
an azeotropic system is unique if ary one of these three 
variables is specified. In order to know the effect of one 
of those variables over the others for the azeotropes, the 
present analysis is divided into two broad classes; 

( 1 ) Prediction of 1 knowing P and x. 

(ii) Prediction of T and x at different values of P. 

The prediction of the P-T-x space for both the classes 
is based on Redlich-Kister' ( 3 . 1 ), Van-Laar ( 3 . 2 ) and Seatchard- 
Hammer (3.3) equations of third order for excess Gibbs free 
energy. The eductions are as follows; 

G E = RTx-jXr, I a 1 + B 1 (x 1 - xp) 1 (3.1) 

f 2 2 

T-, ’ A q x. x 9 B 9 x 9 x. 

g e =rt| — 4 — - — — — + ; ( 3 . 2 ) 

I A v2 f , B 2 v 2 f 

< ^B 2 X 1 ~ x 2^ ‘ r Ag 2^ ’ 

-r? 1 O " p 

G E = RT* Z^(A 3 + 2Z 1 (B 3 -A 3 )) + x 2 Z^ + 

VL ? * 1 

2 Z 2 (A 3 VL" “ B 3^ j 


(3.3) 
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where, 


x. 


x 1 + x 2 


vl 2 

VTAj' 


x 2 VL 2 /VL 1 

vn~ 

X 1 +X 2 vi^ 


For the regular solutions these equations reduce to the 
following forn ; 

G E = A 12 x^ 2 (3.4) 


Here the constant A^ 2 is independent of temperature. In the 
present study the temperature dependence of the constants m 
the equations (3.1) to (3.3) is assumed to he as follows; 


T h A ih 1 ~ ’AfPij 

= Wl 2 


(3.5) 

(3.6) 


where, i = 1, 3. 

This temperature dependence of the constants is based 

on the regular solution theory which essentially presumes the 

E 

product T In . 1 constant. From these expressions of G- one 
can obtain the activity coefficients using the following 
relation ; 


In, 1 =- 


1~ B 
AT / 1 


1_ 

RT 


(n™ Ct E ) 




n. 


n 1 + n 2 


(3.7) 


which can be used m predicting the P-T-x space of the 
azeotropes. The expressions thus obtained tor activity coeffi- 
cients are as given below; 
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For R-K Equation 
In 


y -j = ! Ai +B-j ( -Xg ) * + ^2 ! (og-x-j) + 3^(6x^Xg - 1) * 

L*- * 

(3.8) 

> ] - , -j 

In ^‘2 = x -] x 2 i +B 1 ^ x i “ x 2 ^ X 1 1 ( x 2‘“ x 1 ) + ®-j( ^ x i x 2 “ 1 M 


For Van-Laar Equation s 

2 


In y 


h. 


Xr 


1 fh. % 2 

( BT X 1 + x 2 ^ 


(3.9) 


( 3 . 10 ) 


„ Bp x 1 

In )> 2 = £ 1 


B 2 2 

(x 1 + x 2> 


For Scatchard-Hammer Equation s 


In 


In 


= 4 I ^ 4 2Z 1 (B 3 - A3) 


VL. 


1 


,2 i 


VL 2 h 


V 2 = Zi j b 3 + 2Z 2 (a, b 3 ) j 


(3.11 ) 


(3.12) 

(3.13) 


For estimation of activity coefficients, it is needed 
to evaluate the constants of these equations hy rearranging 
the respective equations. It requires the values of the 
activity coefficients at some known pi 'assure -temperature and 
composition. For this the activity coefficients are calculated 
assuming the vapour phase m equilibrium viath the liquid phase 
to be a nonideal mixture. The equations are; 

m y = m ) + fa j (B^-viptp-iP) + p r 12 1 | j 

(3.14) 
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In / 2 - ( p / p 2 x 2^ + ST' ( B 2 2 -PPl 2 ^ ^ ^ + P ^12 

(3.15) 

where f^ 2 = 23 -i2 ~ 3 n “ B 22 

Here the interaction virial coefficient B., 2 is assumed to he the 
harmonic meanof the individual virial coefficients, and 


1 


B, 


22 ’ 


i.o, 


2_ 

Bh 


1 

B- 


+ 


B, 


( 3 . 16 ) 


1 2 "1 1 "22 
The equations (3.14) and (3.15) require (i) P-T-x data instead 
of P-T-x-y data because x = Y, (li) pure component vapour 
pressures P° , P 2 which are estimated by the Antoine equations, 
(in) the second virial coefficients B^ and B 22 which are 
predicted by employing the Kreglewski's double square well 
potential, and ( iv) the pure liquid molar volumes VL^ and VL 2 
calculated by dividing the molecular weights by respective 
specific gravities because of the assumption that YL^ and VL 2 
do not charge ’"lth temperature end pressure significantly 
and also VT^ and VL 2 m comparison wibh and B 22 are quite 
small. Details of the equations used for the vapour pressure 
and the second virial coefficients evaluation are as follows; 
(i) An to inno equation for vapour pressure calculation; 

log 1Q P = A - B/ t + c (3.17) 

Constants A, B and c are obtained fmm the ' iterature for a 
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particular pure component under consideration* 

( 11 ) Kreglewskis model gives the following equations for 
estimation of the second virial coefficients, and B 22 1 

b li -u At - - -su/kT- 

1F" = 1 " ( R rU(e -D-CEg-RpCe -1 ) 

(3.18) 

where, and are given by the following expressions; 

&l = ! 1 + _s __ | (3.19) 

' (y *) ' ^ 

f V* / , 3 

= » 1 + ( y|) 1/3 j (3.20) 

t 

Here, V* is the liquid molar volume of methane at reduced 
temperature, Tr = 0 . 6 ? s is a constant and is equal to 0.21 
rounded off for all the substances and V* is the liquid molar 
volume of the compound under consideration at Tr = 0.6. The 
potential energy term u/kT is estimated by; 

-u/kT = — E* x 0.239 (3.21) 

where, A is the density and is assumed to be constant at 0 . 32 , 

1 1 

E* is the energy term and can be obtained from the literature 
for a particular compound, 0,239 is a conversion factor used 
and T is the temperature at which the virial coefficients are 
to be calculated. The packing factor, b is estimated by the 
equation ( 3 . 22 ), l.e.s 

b = 0.77 + 0.1325 (T c ) 1//3 


( 3 . 22 ) 
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where, f c is the critical temperature of the pure compound 
under consideration. 

From the knowledge of the activity coefficients at one 
point, the constants m the equations (5.1) to (3.3) can be 
obtained by the following expressions ; 

i i 5 / .2 I— 2 \ — 


y ^ 4 * *■* Xr > ) 


(3*23) 


=i - ( 


In y . 


'2 \ // X 1 X 2 + 5 x 1 x 2 " *2 


x 1 - x^x 2 - 5x^2 


, _ , o *2 l0 S ^2 ,2 

% - log ,, 1 (1 + — xj-lbY 

x 1 log >' 2 

E 2 = log > 2 (1 + J— -i) 2 

x 2 log y 2 


(3.24) 


(3.25) 


(3.26) 


In y 


? ViJ 2 

2z i B 2vir 


B 5 (1 -- 2Z 2 


J 2 VL„ 


(3.27) 


Z 2 (1 - 2Z.j ) 

2z 1 n^vig 

( d- 2ty 


) - ( 


2Z^" z 2 vl 2 /vl 1 

1 - 2Z 2 

2 ? 2 YLg/Vl^ 


(3.28) 


Now for predicting the P-I-x space, C- values are 
calculated using equations (3.1), (3.2) art (3.3). The azeotropic 
temperature can be obtained from the expression; 
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Ct E = RT (x 1 In > ) 1 + x 2 In V g ) (3.29) 

Rearranging equations (3.14) and (3.15) one gets; 

Y^p" 


(Vl rBl1 ) (P-P°) - 4 1 2 , 2 


and 


PY, = p“ X, y, exp! — i— U-_ ^ — — — ;(3.30) 

U ■*"’ 

Using equations (3.30) and (3.31 ) in (3.29) we have: 


PY, 


G E = RT (x- In ~ — - 

1 -rP _ 


PYn 


P^ exp 


+ 


*2 


In 


P 2 X 2 exp U 2 


-) 


(3.32) 


_ „ (^i - Bh)(p-^) - 3 12 4 p 

v/t iGT6y U 1 — r ""' 1 RT ' 1,1 ri "" " ' ' 

, „ (YI 2 - %> (P-I§) -r 12 4 p 

and u 2 ^ 

equation (3.32) may be further simplified to: 

Y Y 

G E = RT (x„ln ~r f xb In ~ + x 1 In ~ + x 2 In ~ 

(3.33) 

(3.34) 


T 


1 


"1 


or 


- x^ - y^v z ) 

g e = ET (V + V - u) 


where, V = x^ln — + x? In 


1 


P_ 

pC 


2 


Y Yr 

W = x^n — + x 2 In 


U = x 1 U 1 + x 2 U 2 

Hence from equation (3.34)? the azeotropic temperature 


is obtained from the simplified expression: 
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T 


g e. 

Tf + W - u) 


(3.35) 


Equation (3.35) is implicit m T. Hence T can be 
obtained by an iterative procedure w th the first approximation 
on the right hand side taken as T = T x 1 + T 2 x 2# Further, 
computer calculations show that Y is very sensitive to tempera- 
ture and even a little off approximation for T on the right 
hand side makes the expression (Y+W-U) negative, yielding 
negative values of temperatures. To obviate this difficulty, 
v is expressed m a different form. The claussius-Clapeyron 
equation can be expressed interms of the entropy of vaporization 
as follows: 

RT in P/l£ = ,6S 1 (T i - T) (3.36) 


In this expression it is assumed that the molar heat of 
vaporization is constant over the temperature range under 
consideration. Here I is the pure component boiling point 
at pressure P and /)_ S- is the entropy of vaporization. 
Substituting this m equation (3.54), the following expression 
results : 

G E = x 1 / S 1 (T.J-T) + x 2 As 2 (T 2 -T) + RTW - RTU (3.37) 

Hence T can bo obtained from equation (3.37), whose explicit 
form is as follows: 


T = TlXl + * 2 * 2 
where, £S 0 = x^ /,\S^ 


/'- S 2 ' C E 
A S o " L % 

+ ^2 / a ^ 2 


RTW 
A S. 


RTU 

A S o 


4 - 


(3.38) 
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This is the final expression which will he used for the 
prod. ict ion of T m the azeotropic P— T— x s^ace. 

for prediction of T knowing P and x, as a first approxi- 
mation, T is obtained as (T.^ + T 2 x 2 ) . The entropies are 
calculated from equation (3.36). The value of G E will be cal- 
culated from equations (3.1) to (3.3) separately. Here TJ is 
the expression for nonideality in the vapor phase. The term 
involving V m equation (3.38) is neglected because of the 
reason that for an azeotrope x = Y. Also m equations (3.14) 
and (3.15) Y^ and Yg are replaced by x^ and x 2 respectively. 

Hor the prediction of T and x at different values of P, 
all the above equations are used. The procedure for this type 
of prediction is as follows! 

Suppose P f is the pressure at which azeotropic compo- 
sition and boiling temperature are desired. Pirst of all 
corresponding to pressure, P' approximate azeotropic boiling 

7 

temperature is obtained by the method given oy Skolnic 
(Chapter 2). lit cose the components belong to the same group, 
then for estimition of T^ and Tp from nomographs, normal boiling 
temperatures of the pure compounds are used instead of the 
normal azeotropic boiling point. Wow using this approximate 
azeotropic temperature, the values of the activity coefficients 
are obtained using equations (3.14) and (3.15) m which the 
last term containing composition of the vapor phase, Y is 
neglected, because of the fact that the term, P 2 Y^ ls very 
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small compared to the other terms. Further, equations (3.8) 
to (3.13) are used to calculate the activity coefficients at 
a number of assumed compositions, x and then they are matched with 
the activity coefficients obtained at the approximate boiling 
temperature, T using equations ( 3. 1 4)-(3 .15) as mentioned earlier. 
Hence corresponding to pressure, P' approximate values of 
temperature, T and composition, x are obtained. Finally these 
initial values of T and x are further improved using the 
iterative procedure described earlier, checking m between whether 
x approximately equals Y and also whether the term W m equation 
(3.38) goes to a very small value. 

For these two type of predictions as stated above 
computer programs have been developed to determine the P-T-x 
space of the azeotropes. 


*•*** 



CHAPTER 4 


RESULTS AND DISCUS SION 

In this chapter the results of two ways of predicting 
azeotropic data namely (i) azeotropic temperature knowing 
azeotropic pressure-composition and (ii) azeotropic temperature 
and composition at different pressures will he presented. Por 
both types of prediction azeotropic (P-T-x) data are necessary 
at one point. Analysis of the methods and all necessary equations 
are presented m Chapter 3. The first kind of prediction is 
essentially to test the prediction methods for their accuracy 
as it is known that both temperature and composition are unique 
at a fixed pressure. The equations used for the two kinds of 
predictions are essentially the seme. The Gauss-Siedel iterative 
technique is employed for solving these equations on computer and 
is briefly described m Appendix A. The computer program for 
the two types of predictions are given m appendix B. 

In order to tost the validity of the prediction methods, 

the systems chosen are diversified in their physical and 

structural properties. They belong to the groups ; aromatics- 

paraffms, aromatics napthenos, hydrocarbon-alcohol, 

e s t © r - alcohol and alcohol-water. Por all these systems 

azeotropic data were available m the literature. The physical 
1 2 

properties needed for calculation of the activity coefficients 
are presented m Table 1. These properties have been predicted 
using the methods mentioned m Chapter 3, wherever experimental 
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data were not available. The constants in the equations of 
Redlich-Kister, Van-laar and Sea tchard -Hammer, have been 
obtained from knowledge of the activity coefficients at known 
compositions and they are presented an Table 2. The extent of 
nonideality and assymetry of the systems are manifested by the 
magnitude of these values, 

Tho azeotropic temperatures predicted using the first 
type of prediction are presented in Table 3. Comparing these 
temperatures with the experimental values it is seen that the 
agreement is very good for nonpolar-nonpolar as well as polar- 
nonpolar systems. Relatively large devietions observed in the 
systems, comprising polar compounds may be attributed to the high 
nonideality caused by the association between the two components 
and strong dipole-dipole interactions. The values of activity 
coefficients, vapor pressures, virial coefficients, and 
constants of the three equations used at these predicted 
azeotropic temperatures are presented m Appendix - A, 

The azeotropic temperatures and compositions obtained 
using the second type of prediction arc presented m Tables 4 
to 6. Por this prediction approximate values of azeotropic 
temperatures, figures (4.7) to (4.12) and ranges of compositions 
are first obtained by the method described m Chapter 3 and 
are presented m Tables 7 and 8 respectively. The computer 
program needed for approximating the azeot- opic temperatures 
is given in Appendix B. Prom the comparison of the predicted 
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VALUES OE THE CONSTATS IN R-K, V-L, S-H EQUATIONS AT AZEOTROPIC 

TEMPERATURE 
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TABLE 3 

COMPARISON OR OBSERVED AND PREDICTED AZEOTROPIC TMIERAlURE 
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FIG .4 7 APPROXIMATION OF AZEOTROPIC TEMPERATURE FOR 
. SYSTEM ETHYL ALCOHOl-ETHYLACETATE 









4.9 APPROXIMATION OF AZEOTROPIC TEMPERATURE FOR SYSTEM 
HEXANE-BENZENE. 



FIG- 4.10 APPROXIMATION OF AZEOTROPIC TEMPERATURE 
FOR SYSTEM METHYL ALCOHOL-BENZENE 






1000/ tc ■+ 230 ► 

FIG 4.12 APPROXIMATION OF AZEOTROPIC TEMPERATURE 
FOR SYSTEM . ETHYL ALCOHOL-WATER . 
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TABLE 7 


APPRO XQ1ATION OR AZEOTROPIU 


TEM PERaTURE S 


SI. 

No. 


System 


Known p 
(P-T-x) mm Eg 




T az ( a P P ) 

°K 


Ethyl alcohol + 
ethyl acetate 


P= 760.0 423.O 

x , =0.4610 220.0 

T az =344.96 774 


337 . 5 ' 

323.0 

302.8 


333*4 333.1' 

316.8 320.1 * 

294.0, 301 .A'' 


2. Benzene + 
Cyclohexane 


3 * Hexane + 
Benzene 


P=759.0 
x. =0.5020 
T a ;=350.56 

P= 760.0 
x. =0.9630 
T az -341.96 


4* Methyl alcohol+ 
Benzene 


P=678.4 
x, =0 .6080 
T az - 3 20.16 


5. 


Ethyl alcohol + 
Methyl cy cl ohexan e 


P=120 .2 


x 1 =0.5400 
T a ’=303.16 


P= 760.0 
x, =0.0940 
T a ;=35i,26 


600.0 

206.0 

345.7 

316.0 

348.2' 

317.9 

349.8 

318.6. 

735 .0 

cb 

• 

O 

KA 

352.1. 

342.2 

366 .0 

320. r 

331.2 

324.8' 

41.3 

292. it 

295.7 

295-3 


380.0 

334 8 

354.8 

335.5 

190.0 

320. CL 

338.4’ 

321.7 

95.0 

306.6 

323.6 

308. C 

50.0 

295.2 

311.3 

297. 81 


6. Ethyl cilcohol + 
Water 
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azeotropic temperatures and compositions (Tables 4 to 6) it 
is observed that the agreement is very g^-od for all systems 
excepting the polar-polar mixtures. 

The equilibrium liquid and vapor compositions at various 
predicted azeotropic temperatures have been calculated for all 
systems and arc presented m Table 9 • These equilibrium data 
are plotted m Figure (4.1) to (4.6). 

From the analysis of the results it is observed that 
all the three correlations used for calculation of the excess 
free energy yield equally good agreement of the predicted 
azeotropic data with the experimental values, although the 
inverse temperature dependence of the constants m these 
correlations may not be valid for highly polar systems. 

The first approximation of the azeotropic temperature 
should be done carefully. This is because equation (3.38) is 
implicit in T and is highly non-linear so it will have more 
than one solution. Hence a little off approximation outside 
a range will result false convergence. 

In the Gauss-Scidel iterative method , the convergence 
is relatively slow. Although fast converging methods such 
as the Newton-Raphson method are available, all these methods 
require analytical derivatives of the function used. In the 
present study the function is quite complicated and hence it is 
a tedius task to evaluate analytical derivatives of this 


function. 
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TABLE q 

equilibrium data at IIPEERMT pressures 


Si. P 

No. System mm jj g 


1. Ethylalcohol + 423.0 

etbylacetate 


220.0 


77.4 


2. Benzene + 
cyclohexane 


600.0 


206 .0 


3. Hexane + Benzene 735*0 


4* Methyl alochol + 366 .0 

Benzene 


V»° K 

CL 4J 

*1 

Y 1 

329. 9< 

0.3000 

O.3335 


0.3500 

0.3681 


0.4000 

0.3993 


0.4500 

0.4280 


0.5000 

0.4548 

314.6' 

0.2500 

0.2758 


0.3000 

0.3112 


0.3500 

0.3425 


0.4000 

0.37C7 


0.4500 

0.3964 

293.K 

0.1500 

0.1700 


0.2000 

0.2106 


0.2500 

0.2454 


0.3000 

0.2756 


0.3500 

0.3021 

343.2 

0.4000 

0.4243 


0.4500 

0.4624 


0.5000 

0.5004 


0.5500 

0.5386 


0.6000 

0.5773 

313.7 

0.4000 

0.4236 


0.4500 

0.4603 


0.5000 

0.4969 


0.5500 

0.5337 


0.6000 

0.5712 

340.8 

0.8500 

0.9150 


0.9000 

0.9350 


0.9500 

0.9550 


0.9700 

0.9650 


1 .0000 

1 .0008 

313.1 

0.5000 

0.5773 


0.5500 

0.5795 


0.6000 

0,5834 


0.6500 

0.5896 


0.7000 

0.5984 
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Table 7 (contd) 


S'. Syste ” 

P 

mm Hr 

T , °K 
az 

X 1 

Y 1 

5. Ethyl alcohol + 

41.3 

283<3 

0 .4000 

0.5180 

Methyl cyclohexane 



0.4500 

0 .5046 




0.5000 

0.4935 




0.5500 

0.4852 




0.6000 

0.4804 

6. Ethyl alcohol + 

380.0 

334-7 

0 .8000 

0.8181 

Water 



0.8500 

0.8569 




0 .9000 

0.8988 




0.9500 

0.9438 




1 .0000 

0.9919 


190.0 

319.8 

0.8000 

0.8171 




0.8500 

0.8554 




0.9000 

0.8968 




0.9500 

0.9415 




1 .0000 

0.9894 


95.0 

306. 3 r 

0.8000 

0.81 39 




0.8500 

0.8514 




0.9000 

0.8922 ‘ 




0.9500 

0.9364 




1 .0000 

0.98A0 


50.0 

295 -0 

0.8000 

0.8143 




0.8500 

0.8513 




0.9000 

0.8918 




0.9500 

0.9358 




1 .0000 

0.9832 
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1 . 
2 . 

3. 

4. 

5* 

6 . 

7. 

8 . 

9. 

10 . 

11 . 

12 . 

13. 

1 4 . 

15. 

16. 

17. 

18. 
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APPENDIX A 


1. Gauss-Seidel Iterative Method -. 

Gauss- Seidel iterative scheme is used when the equation 
is of the form: 

x - f(x) 

In the beginning of the iterative scheme initial guess 
of the value x is necessary. Next in the process of iterations, 
the newly computed values of x are taken as the next better 
guess to continue the process, Finally we get the solution when 
the difference between the newly computed value and the preceding 
one is negligible. 


2. Vapor Pressures, Activity Coefficients, Virial Coefficients 
and Constants of the ^hree Equations at Azeotropic Temperature 
Composition ? 

Units of P, £}, Pg - mm of Hg 
Units of - °K 

cI/j 

Units of B^, B 22 ~ cc/mole 


a. System ; Ethyl Alcohol + Ethyl Acetate 


P = 425.0 

T az =529.9 

X 1 az 

P° = 505.55 

^ = 1.5742 

B 11 

Pg = 570.50 

„ 2 = 1.1562 

B 22 

A 1 = 0.8555 

A 2 = 0.8245 

H 

B 1 = 0.0520 

B 2 = 0.8889 



0.5970 

- 2527.5908 

- 1807.5570 
0.8041 
0.8048 
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P = 220.0 

T az =314.6 

X 1 az = °*3310 

p ° 

= H4.79 

3 1 = 1.5013 

B^ = - 3025.6215 

*2 

= 200.17 

: 2 = 1.0965 

b 22 = - 2097.4711 

A^ 

= 0.8970 

^2 = 0.8645 

Aj = 0.8431 

B„ 

1 

= 0.0335 

B 2 = 0.9321 

B 3 = 0.8439 


P = 77.4 

T az = 2 93.1 

X 1 az = °*2390 


= 43.83 

*;■, = 1.7534 

B u = -4047.5522 

tP 

*2 

= 73.62 

> 2 = 1.0505 

B 22 = -2658.5524 

A 1 

= 0.9629 

A 2 = 0.9280 

a 3 = 0.9050 

B i 

= 0.0360 

B 2 = 1.0005 

B 3 = 0.9058 

System 

t Benzene + Cyclohexane 



P = 600.0 

T az = 343.2 

x 1 az.= °’ 4 ?80 


= 551.54 

'/•] = 1.0854 

B 11 = - 1041.8023 

n 

= 508.87 

y 2 = 1.1719 

B ?2 = “ 1224.4343 


= 0.4810 

ig = 0.7134 

A 3 = 0.6068 


= -0.1563 

Bg = 0.3629 

B 3 = 0.2841 


P = 206.0 

T az = 313.7 

x 1 az = 0,485 

tP 

■m 

= 186.75 

= 1.1018 

B^ = -1325.0588 

*§ 

= 173.88 

> 2 = 1.1811 

B 22 = -1573.4160 

*i 

= 0.5263 

Ag = 0.7804 

A 3 = 0.6639 

b i 

=-0.1709 

B 2 = 0.3970 

b 3 = 0.3109 

System 

: Hexane + 

Benzene 



P = 735.0 

T az = 340.9 

X 1 az , = °-9 610 

•pO 

r 1 

= 737.24 

y 1 = 0.9971 

B^ = -1339.8401 

pO 

*2 

= 511.50 

) 2 = 1.4211 

B 22 = -1059.8589 
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JL, = 1.2924 

^2 = -0.0435 


= 

4.8939 

B 1 =-1.0867 

B 2 = 0.2210 

®3 

= 

0.2019 

System; He thyl alcohol + Benzene 




I = 366.0 

T az =313.1 


x 1 az " 0 4 5800 

P° = 265.26 

y 1 = 1.3705 

B 11 

= - 

-1597.3142 

P^ = 182.49 

)? 2 = 1.9730 

B 22 

= - 

-1331.6952 

it, = 1.9401 

^2 - 2.0634 

A 3 

= 

1.7668 

B 1 =-0.1127 

B 2 = 1.8313 


= 

1.3029 

System: Ethylalcohol + Methyl Cyclohexane 



P = 41.3 

Az = 28 ^ 


x 1 as S °- 4935 

P° = 23.65 

= 1.7382 

B ii 

= 

-4704.6868 

Pg = 21.59 

)> 2 = 1.9059 

B 22 

= 

-2854.3253 

A 1 = 2.3952 

^2 = 2.6756 


= 

2.3053 

B 1 =-0.2467 

B 22 = 2.1683 

=3 

= 

1.1315 

System: Ethylalcohol + Water 




P =380.0 

T az = 334.7 

x 1 az = °- 8800 

P° = 376.06 

y., = 1.0101 

B 1 1 

= 

-2398.6809 

ig = 160.24 

y 2 = 2.3409 

B 22 

= 

-1115.1555 

A 1 = 1.2022 

^2 = 1.5880 

■S 


0.0700 

B 1 =-0.1999 

b 2 = 1.0040 

B 3 

= 

1.0027 
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P = 190.0 


T 

pp 

Jm 



az 

= 187.73 

V 1 

= 1.0119 


= 78.08 

% 

2 

- 2.4137 

A 1 

- 1.2582 

Ag 

= 1.6619 


=-0.2092 

b 2 

= 1.0507 


319 - 8 *ia z = °-S755 

B-n = -2840.3379 
B 22 = -1288.5237 
^ = 0.0733 
= 1.0494 


P = 95.0 
3^ = 93.69 
rg = 38.0? 
Al - 1.3133 
B 1 - -0.2 184 

P = 50.0 
3^ = 49.22 
?2 = 19.52 
A 1 = 1.3677 
B 1 =-0.2267 


T az = 506.3 
>■-, = 1.0139 
*y 2 = 2.4844 
^2 = 1.7347 
b 2 = 1.0968 

Az = 295.0 
> 1 = 1.0157 
)? 2 = 2.5522 
Aj = 1.8013 
B 2 = 1.1389 


x 1az = 0.8670 
= -3363.7714 
B 22 = -1489.6359 
Aj = 0.0765 
b 3 = 1.0954 

x 1az = °- 8610 
= -3934.7287 

B 22 = -1704.7313 

= 0.0794 

B 3 = 1.1374 



COMPUTER PROGRAM 
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.1ST CF IMPORTANT VARIABLES USED IN COMPUTER PROGRAMS ONE T ’-'ROUGH T u P. r E 

}?# 3{c^3^ %%% sfc^sjc :Je * # sjcsfc ^jjc^^jjssfcjjcsfc $ 2 Qe$:$ sjc^ jje 


VL1,VL2 

RHO,RH.OC f S, VO STAR 

A1,B1 ,C1 

A2,B2 ,C2 

VSTAR 

ESTAR 

TC 

TO, XI ,X2,P 
T1*T2 
VC(1, 1) 

VCC2, 2) 

VC(1, 2) 

GAMA1 
GAMA2 
AT (M) 

GE 

VP1.VP2 

B 

RKA,RKB 
VLA,VLB 
SHA » S FB 


MOLAR VOLUME*” OF PURE COMPONENTS ( L IQU T D 1 ONE ( 'N n ’ r WO . 

KREGLESKI S r ONS T ANTS 

ANTOINE CONS T ANT S FOR COMPONENT ONE 

ANTOINE CONSENTS FOR COMPONENTS TWO 

LIQUID MOLAR VOLUME AT REDUCED TEMPER AT’ ’R E T/ T C=0.6 

ENERGY OF VALORISATION OF COMPONENT A T n FDUCc n T?r »'P. 

CRITICAL TEM°ERA T URE 

OBSERVED A Z2 n TRO°I C TEMP COMPOS IT ION A*'P PRE c S’”>r- 
BQILING POINTS OE PURE COMPONENTS IN W I VT URE T N * 
VIRIAL C0EFF T C IE M T OF PURE COMPONENT n N r 
VIRIAL COEFFICIENT OF PURE COMPONENT T W n 
INTERACTION VIRIAL COEFFICIENT OF THE M T XTURE 
ACTIVITY COE c FI C T ENT OF COMPONENT ONE 
ACTIVITY COE r F IC TENT OF COMPONENT TWO 
APPROXIMATE T EMP r RATURE EVALUA T ED 
EXCESS GIBBS FRE r ENERGY 

VAPOR PRESSURES n F COMPONENTS ONE AND TWO 

PACKING FACTOR 

REDLICH-KIST C R CONSTANTS 

VAN-LAAR CQN*”T AN T S 

SCAT CHARD-HA M MER CONSTANTS 




G I 


^ & $ 4 $ # & ^ & 

£ # 

* PROGRAM t * 

* & 
4% 


PREDICTION OF AZEOTROPIC BOILING TEMPERATURES, US I NG REDL ICH-K 7 STER, 
VAM-LAAR A M D SCATCHARD-HAM M ER EQUATIONS. 


C LIQUID MOLAR VOLUMES ARE ASSUMED TO *E INDEPENDENT OF TEMP. AND P°£ . 
C FOR CALCULATION OF VIRIAL COEFF I r IEN T S KREGLESKI S METHOD I c US r U 

C FOR VAPOR PRESSURE CALCULATION A*'TOI M E EQUATION IS USED 

C UNITS OF R (MM OF HG) (C.C . ) / (GM-M^LE I ( K) 


DIMENSION UUC),VC( 10 ,! 0 ) ,VSTA°( 5 ) ,ESTAR( 5 ) ,TC( 5 ),AT( 10 O).UU( 2 u) 

VL 1 = 46 • 07 / 0 . 7 P 9 
VL 2 = 88 . 01 / 0 . 9 '! 

R= 80 . 06 * 760 . 

1 FORMAT!/, 15 X,*R-K CONSTANTS*, / /, 6 X , *G AMA 1 * , 7 X , *GA MA 2 *, 9 V , *RK A*,’ O v 

i,*rkb*,iox,*t:*,uy,*t2*,/) 

2000 F 0 RMAT(/» 2 F 14 . 4 , 3 F 16 . 4 ,F 15 . 4 ) 

2 FOR MAT ( /» 15 X* *VAN-LAAR CONSTANTS* , //, 6 X, *GAM 1 A 1 *, 7 X ,*GAM a 2 * , 9 X, *VL * 
l*,l CX»*VLB*, 10 X,*T 1 *, 11 X,*T 2 *, /) 

103 FOR MAT (/» 15 X,*SCATCHAR D-HA MMER C 0 NCTANTS*,//, 6 X,*GAM!Al*, 7 v ,*GAM* 2 * 
1 , 9 X ,* 5 HA*, 10 X,*SHB*, 10 X,*T 2 *,i 1 X,*T 2 *,/) 

3 FOR MAT ( 1 H 1 ,//, 10 X,*C 0 MPARI SION OF HBSERVED AND PREDICTED T EMPER*T" 
IRES BY R-K,VAN-LAAR AND SC ATCH ARD HAMMER EQUATIONS*, // ) 

5 FOR MAT (/»I 5 X,*F 0 R ^THANOL* , / , 8 V , *VP 1 *, 9 X, *VC ( 1 , 1 5 * , 7 X , *^ 6 ’ TA*, 6 V ,* 
1 VC( 1 , 2 )*,/) 

8 FOR MAT ( /, 15 X, *FOR BENZENE* »/, 8 Y ,*V D l*» 9 X»*VC{l»l)*»?X,* n E* Y A*, 6 Y ,* 
1 VC( 1 , 2 )*,/) 

320 FOR MAT ( /, 15 X,*F 0 R HEXANE* ,/, 8 X , *VPi *, 9 X, *VC ( 1 , 1 ) *, 7 X , *D C L^A*, 8 X , *” 
1 C ( 1 , 2 )*,/) 

321 FOR M 1 AT(/» 15 X»*F 0 R METHANOL *, / , <> X,*VP 1 *, 9 X,*VC( 1,1 ) *, 7 X,*D’"LTA*,°X, 
1 *VC ( 1 , 2 )*,/ ) 

322 F 0 RMAT(/» 15 X»*F 0 R METHYLCYCLOP c NTAUE*, /, 8 X , *VP 1 *, 9 X , *VC ( 1 , 1 ) * , 7 Y , * 
1 DELTA*, 8 X,*VC ( 1 , 2 )*,/) 

6 F 0 RMAT(/» 15 X» *FOR ETHYL AC ETAT r * , / , 8 X , *VP 2 * , 9 X ,*VC ( 2 , 2 ) * , '’X, *Dc’ T « 
1 *, 8 X,*VC( 1 , 2 )*,/) 

11 FORMAT! / , 15 X,*FOR CYCLOHEXANE *,/, 8 Y,*VP 2 *, 9 X,*VC( 2 , 2 )*," T X,*DELTA*, 
18 X, *VC ( 1 , 2 ) *, / ) 

323 FOR MAT ( / , 1 5 X , *FOR METH YLCYCLOH r XAN c *» / , 8 X,*VP 2 *, 9 X,*VC('*,'’)*, 7 X,* n 
1 ELTA*, 8 X,*VC( 1 , 2 )*,/) 




ICO > FORPAT(/,15X,*FOR R eNZENF*,/ T 3 Y y*VP2*,9X,*VC(;?,2)*,7X,* n e , TA* t o v ,* 

1 VC( 1,2)*,/) 

7 F0RPAT(4F16.1,/) 

22 FOR PAT ( 1H1, / , 50X, SYSTEM- ETHYL ALCOHAL AND ETHYL ACF.TA" r E*» / f 5v v i‘ ) 
1 9 { 1 F* ) ,//,5X»*KNOWN POINT- * , 2X,*P=*,F7.1 ,2X,*TO=*, F6 . 1, R X- *X1=*, F^ 
2.4, 2X,*X2=*»F7.4,/} 

33 F0RPAT(1H1,/,50X,*SYSTEM- BENZENE AND CYCLONE XAN*=* ,/, 5. v , ~ H 1 H* 3 , ' 
1 / , 5 X , * KNOWN POINT— 4 , 2X,*P=*»F7 .1 ,2 Y ,*TG=*, F6. 1 , 2 Y , *X I=*. F” . 4, 2 X , * Y 

2 2=* , F7 . 4, / ) 

314 FOR PAT ( 1H1,/ ,50X,*SYST£M- HEXANE A M D BENZENE* ,/, FOX, 26 PH* ),//, r X , 

1 *KN OWN P0+NT-*,2X,*P=*,F7. 1 , 2 X ,*TO=* , F 6. 1 , 2X, *X1 =* , F7 .4 , 2 V , *X2=*, ^ 
27.4,/) 

316 F0*PAT(1H1,/»50X,*SY2TEM- METHANOL AND BENZ ENE*, /, SOX , 2 R H* ),//, c 

IX,* KNOWN PO+NT-*, 2X, *P =* , F7. 1 , R X , * T 0=* , F6. 1 ,2X , *X1 =* , F7.4,7X, */'»=* 

2 , F7 .4, / ) 

318 FOR PAT( 1H1, / ,50X,* C YSTEM- HETH V LC Y r LDP ENTANE AND B6NZEN c *,/,5 X,2 ° 
1 1 1H *)»//» 5X , * KNOWN POINT— *,2X,*P = *,F7.1,2X,*T0=*,F6.1,2 V , * Y 1 = *, c 7. 
24,2X,*X2=*,F7.4,/) 

319 FORPAT(1HI,/,?OX,*$YSTEM- ETHANOL AND METHYLCYCLOHEX ANE*» / , 50X ,“»7 f 
11H* ),//,5X,*K?!OWN POINT-*, 2X, *^ = * , P7. 1 ,2X, *TO=* , F6 . 1 , 2X , * V I =* , F” . 

2 , 2X ,*X2=*, F7.4 , / ) 

77 FOR PAT(/,5X,F7.4,5 y ,F7.4,5X,F9.4,5X,F8.4,5X,F8.2,5X,FG.'', /} 

99 FORPAT( F4.1,F5.3,F C .3, F6.1 ) 

2 22 FOR P AT ( F3 .1 , F4 . 2, F4. 2 , F4- 2 ,F5 , F6 . 2, F 6. 2 , F6 . 2 ) 

109 FOR PAT ( F7. 5 » F£ . 3 » F 7 .3» F7.5,F8.”’»F7 .3) 

666 FOR PAT ( F5.3,F 7 .l,F6.2) 

66 FOR PAT ( /,5X,F3.2,5 Y ,F9.2,5X,F8.2,5V,F9.2,5X,F9.4,/) 

4 FOR PAT ( / »5X»* CHECK P0INT*,SX,*°=*,F7.1,2X»*T0=*,F6.1,2X,* V 1=*,F - ’. * 
l,2X,*X2=*,F7.4,/3 

104 FOR PAT{/,2X,*APPR0XIMATE TEMPE°ATU®E E VALUATED=*, F8.2 , / ) 

88 FORPATl /,9X,*C1BS T^HP* , 10X ,*RK T E * , 1 2X , *VLTE*, 12X , *SHTE* , / 3 
D06 OOMP = 1 , 6 
M=1 
MN= C 
I N=0 
MA=C 
PPI=1 

REAC222»F» RHO , RHOC , S,VOSTAR,RK' r , VLT , SHT 

R6AC109,A1,B1,C1,A?,B2,C2 

D03 3 3 11 = 1,2 

REA 0666, VST AR( II ) , R STAR{II },TCUI) 

333 CONTINUE 

I F( PP.EQ.1)MQ=6 
I F( PP. EQ.2 ) MQ=4 
I F( PP.EQ.3 )MQ=2 
I F( PP. EQ.4 ) MQ=2 
I F( PP . EQ . 5 ) MQ = 3 
I Ff PP. EQ. 6) MQ=2 
D02 C0PI=1 , MQ 
DIF P IN=100. 
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jj=o 

NN=0 

REA C99 » TO * X 1 *X2 ,P 
UU( X) =T0+273 . 16 

T1=B1*( 1 ./{ A]-ALOG(P) /2-303) )- r l+2' 7 3.16 
T2=B2*(1./(A2-ALGG(P)/2.3 j 3) ) -02+2^3. 16 
U(H ) =TO+273 .16 
I F( K-l )10,j Ot’OO 
IQ D0400 1=1,2 

CAL L VIRIAL (U,RHO,$,VOSTAR,VST , 'R,E < :TAR,VC,TC,M,I) 
4 CO CONTINUE 

DELTA=2.*VC(1,2)-VC(1, 1)-VC(2,'’) 

T=TC 

CALL VPC ( T» A 1 ,B1 ,CT , A2 , B2 , C2 , V°1 , V D 2 ) 

I Ft XP. EO.l ) PRINT22, P,TO,Xl ,X2 
IF( XP . EQ. I ) GO 1044 
101 VL1 =78. /O. 879 
VL2 =84. /G.779 

I F{ XP.EG.2)PRINT33,P»TC,X1 ,X2 
I Ft XP. EQ. 2 ) GOT044 
VL1 =86. /O. 659^7 
VL2=78./0. 879 

I F( XP.EQ.3)PRINT31*,P,T0»X1,X2 
I F( XP. EQ. 3 ) G0T044 
VL1 =32. /Q.7924 
VL2=78./0.879 

I Ft XP.EC.4)PRINT316,P,T0,X1,X2 
I Ft XP . EQ.4 ) G0T044 
VLl =84. /Q.750 
VL2 =78. /O. 879 

I Ft XP.EG.5) PRINTS 18, P, T0,X1, X2 
I Ft XP.EQ.5) GCT044 
VLl =46.07/0.789 
VL2 =98. /O. 769 

I Ft XP.EQ.6) PRINT31°,P, TO, XI, X2 
44 AA=AL0G(P/VPI) 

BB= F/ ( R#U ( M ) ) 

CC= (VC(1,1)-VL1)*( D -VP1) 

DD=P*DELTA*X2*X2 
I Ft XP. EQ. 1 1 PRINTS 
I Ft XP.EQ. 2) PRINTS 
I Ft XP.EQ.3 ) PR1NT320 
I Ft XP.EQ.4) PRINTS 27 
I Ft XP.EQ. 5)PRINT32? 

I Ft XP.EQ.6) PRINTS 

PRI NT 66 »VPI»VC(1,1) , DELTA, VC (1,2), VLl 
AAA=AL0G(P/VP2) 

BBB=F/(R*U( M) ) 

CCC = (VC(2,2 )-VL2)*(P-VP2) 

DOD = P*DELT A*X 1*X1 



IF ( PP.EG.i) PRINT6 
I F{ NP.E0.2) PRINT11 
I F{ f'P.EQ.3)PRINT1000 
IF( NP.EQ.4)PRINT100o 
I F ( NP. EG. 5 ) PR INTI': "0 
I F( ^P.EQ.6)PRINT10O0 

PRI M66, VP2,VC(2,2) , DELTA , VC ( 1 , 2 ) , VL2 

GAM Al-E XP( AL0G(P/VP1)+ (F/ ! R*U(M) ) )*( (VC( 1,1 ) -VL1 )*{ P-V^l ' +P*DLL T 
1 *X2 *X2 ) ) 

GAMA2=£XP( AL0G(P/VP2)+ (F/ (R*U ( M ) ) )*< (VC(2>2) -VL? > *( P-V n 2 ’ +P*Dc< T 
I*X1 *X1 ) ) 

CAL L RK{ GAM A 1, GAM A?, XI ,X2, FN, R* A , R*B , GE, R, AT , M ) 

PRI NTL 

PRI KT7T » GAMA I , GAM A? f RK A, RK8, Ti ,T2 

CAL L VANLARI IN,GAM« 1, GAMA2 ,X1 ,V2,VLA, VLB, GE ,R * AT , M ) 

PRI NT2 

PRI NT77 ,GAMA1 ,G AMA* > , VL A , VL B ,TI , T2 

CALL SCAHAM (GAM A] , GA MA2 , XI , X’ , MA , SH A ,SHB, GE, R, AT ,M , VL 1 , V 'L2 ) 

PRI NT1P3 

PRI NT77,GAHA1,GAMA' , ,SHA,SHB,T1,T2 
N=N + 1 
PRI NT3 
GOTO 200 

100 AT( F)=Tl*Xi+T2*X? 

PRI NT4,P,TQ,X:,X2 
PRI ATI 04, AT { M ) 

PRI NT88 
AT( M-l )=RKT 

555 RKA =RKA*AT( H-' )/AT(M) 

RK8=RKB*AT(M-L)/AT(M) 

AT( P-1)=AT(M) 

R KT = A T ( M ) 

ZZI =RKA 
Z Z2 =RKB 

CALL RK(GAMAliGAMA?,Xl»X2, MN, R*A , RKB,GE , R, AT , M) 

444 U(M)=AT(M) 

D05C0 1=1,2 

CALL VI RIAL (U , RHO , S , VO STAR , VST «R , E ETAR , VC, TC , M, II 
500 CONTINUE 

DELTA=2.*VC U *2 J-VC ( 1, 1 J-VCC 2 .*>) 

T=A T ( M ) — 273 .16 

CALL VPC(T, A1,B1,C1, A2,B2,C2, VP1 ,VP2) 

V 1= ( ( VL1-VC (1,1) ) *(P-VP1 )-DELT«*X2*X2*P)/ (R*AT(M) ) 

V 2= ( (VL2-VC (2,2) ) * ( P-VP2 ) -DELT »*X1*X1*P ) / ( R*AT( M) ) 

V=X 1*V1+X2*V2 

DELS1=R*AT(M)*AL0G(P/VP1)/(T1-AT(M) ) 

DEL S2=R*AT ( M ) *ALOG f P / VP2 ) / (T2-*T ( M ) ) 

DEL SO=Xl*DELSL+X2* n ELS2 

AT1=T1*X1*DELS1/DELS0 

AT2=T2*X2*DELS2/DELS0 
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AT3=GE/DELSQ 
AT4=R*AT ( M ) *V/D£LSO 
AT ( F)=AT1+AT2-AT3-AT4 
BT= AT I M ) -2 7 3 » 16 
DIF F = ABS ( BT-TO ) 

I Ft CIFF.LT.DIFMIN ) GOTO 111 
10C4 I F( NN-1 1300,700,1001 
111 E£= AT ( M ) 

0P=VCU,1) 

00= VC (2, 2) 

DR=VC( 1,2) 

0S= DELTA 
DIF F IN = D IFF 
I F( JJ— 1 ) 555 » 55 , 102 
300 RKT E=EE 
OT=CP 
, 0U= CQ 
OV=CR 
OX=OS 

ATI F)=T1*X1+T2*X? 

ATI F— 1 ) =VLT 
DIF F IN= 100. 

JJ= JJ+1 
F!M= FM+1 
NN= NN+ 1 

55 VLA = VL A*AT I M- 1 ) / AT ( M ) 

VLB = VLB* AT (M-1)/AT(M) 

ATI F— 1 ) =AT ( M ) 

VLT =AT IF) 

III =VL A 
ZZ2=VL8 

CAL L VANLARI I N, GAM A 1 , G AMA 2 , XI , v 2 , VLA , VLB ,GE t R * AT » M ) 

GOT C 444 
7 CO VLT E=EE 
PT= CP 
PU= CQ 
PV=CR 
PX=CS 

ATI F)=T1*X1+T?*X2 
ATI F-l )=SHT 
DIF F I N=1 00. 

JJ= JJ + 1 

NN=hN+ 1 

102 SHA = SHA*ATI M-D/ATIM) 

SHR =SHB*AT I M-l ) /AT ( M) 

ATI F-l )=AT( M ) 

SHT =AT I F! } 

SHT E = EE 
QT= CP 
QU=CQ 



QV= CR 
QX= CS 
ZZ2 =SHA 
ZZ2 =SHB 

CALL SCAHAM f G AMA 1 , GA MA 2 , X 1 , X* 1 , MA, SHA , SHB , GE , R, AT , M, VL% 'T2 ) 
GQTC444 

1002 PRINT?, UUIM) , RKTE ,VLTE,SHTE 

PRI NT2000,RKA,RKB,nT,OU,OV»OX 
PRI NT2000, VL A , VLB , PT, PU, PV ,PX 
PRI NT2000»SHA,SHP »OT»QU»QV*QX 
2 CO CONTINUE 
6C0 CONTINUE 
PRI NT9 

9 FOR NAT ( 1H1 5 
STOP 

END 

C THIS SUBROUTINE CALCULATES THE VAPOUR PRESSURES OF THE CO M P n NENTS 
SUBROUTINE VPC ( T , AT , B 1 ,Ci , A2 , , C2 , VP1 , VP2 ) 

VP1 =EXP(2.303*(A1- R 1/(T+C1 ) ) ) 

VP2=EXP<2.303*( A2-B2/I T+C2)) ) 

RETLRN 

END 

C THIS SUBROUTINE CALCULATES THE V^RIAI. COEFFICIENTS 

SUBROUTINE VIRIAL (U,RHO,S, VQST ft R,VSTAR,ESTAR,VC,TC»M»I ) 

DIMENSION U ( 1 ),VCI10, 10) , VSTA P I5) ,E ST AR ( 5 ) , TC ( 5 1 
J=2 

B=0. 77+0.1325*1 TCI T ) )**( 0.334) 

Rl= (l.+S/I VSTAR (I ) ) ) **3 
R2= (l.+VCSTAR/IVSTARI I )) )**3 
R=1 .987 

C UNITS OF R - G.CAL,C.C./GM.-MO! E,QK 

RP1 =( RHO*ESTAR I I ) ** . 23 9 ) / I R*U I M ) ) 

RP2=S*RP1 

WST AR=VSTAR{ I 3**3 

VCI I, I ) = B*WSTAR*( 1.-IR1-1. )*(F v P(R°l)-l.}-( R2-R1 )*IEXPI°P' , }-1.)) 
IF ( I. LT. 2 ) GOTO 9 99 

VCI I-1,J) = (2.*VC(I-1,I-I)*VC( 1,1) )/(VC(I,I )+VC{ 1-1, I — 1 ) 1 
999 RETLRN 

END 

C THIS SUBROUTINE IS FOR R EDL I CH-K T STE P EQUATION 

SUBROUTINE RK I GAMA'* ,GANA 2 , XI , X** ,MN ,RKA ,RKB , GE , R, AT , M ) 

DIMENSION AT 1 1 00 ) 

MN= NN+1 

IF! MN.GT . 1 ) G0T050 
B1=X2*X2 

Cl= IX1*X1*X2+ P .*X1*X2*X2-X2) 

B2=X1*X1 

C2= (Xl-Xl*X2*X2-5.*Xl*Xi*X25 

# RKB = I (ALOGI GAMA1)/R1)-(AL0G(GAVA2)/B2) )/ I Cl/B 1-C2/B2 ) 

RKA=( AL0GIGAMA1 )/Bl )-Cl*RKB/Bl 


66 



67 


GOTC888 

50 GE=R*AT( M)*(RKA*X:*X2+RKB*X1*X' , *(X'?-X2) ) 

888 RETLRN 
END 

C THIS SUBROUTINE IS FOR VAN-LAAR EQUATION 

SUBROUTINE VAf JL AR ( T N , GAMAI , GAM *2 , X 1 , X2 , VLA , VL H ,G C , R, AT, »» ) 

DIMENSION ATCOO) 

I N= IN+I 

I F( IN.GT.DG0T060 

VLA = (AL0G(GAMA1) )*(!. + (X2*ALQGf GAM 6 2) ) /{XI *ALOG(GA MAI ) > )**2 
VLB = ( ALOG{ GAMA2 ) ) *{ 1.+ ( X1*AL0G fGAMfi 1 ) ) /{ X2*AL0G ( 0AMA2 ) ) 1**2 
G0TC777 

60 GE=R*AT(M)*{X'+VLB*X2/VLA)*{XH(Xl+X2*VLB/VLA) )*{ ( VLB*X' , / ,, LA ) / { Y 1 + 
1X2* VLB/ VLA) ) * VL A 
777 RETLRN 
END 

C THIS SUBROUTINE IS FOR SCATCHARD-HAM V ER EQUATION 

SUBROUTINE SC6HAM{^AMAi»GAMA2»Xl»X” , »MA»SHA T SHB »G C »R»AT»*% v t. 1 » VL"’ ) 
DIMENSION AT { 100 ) 

Z1=X1/{X1+X2=MVL2/VL1) ) 

Z2=X2*(VL2/VL1)/(X1+X2*(VL2/VL'’ ) ) 

MA= MA + 1 

IF( MA.GT.DGOTO70 

SHB = (AL0G(GAMA1 )/( 72*7.2* (1 .-2.*Z1) )-ALOG ( GAMA2 ) / 1 2 . *Z1* 7 1*Z 2*VL' , / U 
1 LI) )/( (2.*Zl*VLl/VL2/( i.-2.*ZD )-( (1.-2. *Z2 ) / { 2.*Z2*VL2/V’ 1) ) ) 

SHA =ALOG( GAMA2 ) / { 2. *Z1 *Z1*Z2*VI 2/ VL1 )- SHB* ( ( 1.-2 . *Z2 )/('’. *Z2*VLV W 
1 LI) ) 

GOT Cl 05 

70 GE=R*AT{M)*{XI*(Z2*Z2*(SHA+2.* 7 1*(SHB*{VL1/VL2) — SNA) ) )+ v 2*f Z1*Z 7 *6 

1 SHB + 2. *Z2* ( SHA* ( VL? /VL 1 )- SHB ) H ) 

105 RETLRN 
END 

CENTRY 
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^PROGR^M I T * 
# * 
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PROGRAM FOR OBTAINING APPROXIMATE °ANGR OF AZEOTROPIC COMPOSITION BV 
MATCHING ACTIVITY COEFFICIENTS AND PRODUCT OF ACTIVITY COEFFICIENTS M*'D 

VAPOR PR C SSU°ES 


C LIQUID MOLAR VOLUMES ARE ASSUMED TO B E INDEPENDENT OF TEM n . AND PRE . 

C FOR CALCULATION OF VTRIAL COEFF I r IEM T S KREGLESKI S METHOD I ** USED 

C FOR VAPOR PRESSURE CALCULATION ANTOINE EQUATION IS USED 

C UNITS OF R (MM OF HG) (C.C .) / (*M-MnLEM K) 


DIMENSION U ( TO ) , V C TAR ( 5 ) , EST A° ( 5 ) , TC ( 5 ) , VC ( 10 , 10 ) , AT ( 1*0 ’ 

1 FORMAT! 1H1,/, 50X**SYSTEM- ETHYL AL r OHAL AND ETHYL ACETATE*,/, 5 w v *‘ 7 
19(1 H* ) ,//t5X»*KN0WM P0INT-*,2X,*P=*,F7.1,2X,*TQ=*,F6.1,''X.*X1=*, F’ 
2.4, 2X, *X2 = * , F7 .4 , / ) 

2 FORMAT (1H1,/»50X»*SYSTEM- BENZENE AND CYCLOHEXAN c *, / , 50V , N ( IH*) , / 
1/,5X,*KN0WN POINT-*, 2 X ,*P = *, F7 . 1 , 2 V, *TG=*, F6. 1 , 2X , *X 1=*, F^.4, 2X,* V 
22=* ,F7 .4, / ) 

3 FOR MAT ( 1H1, / , EOX , *S Y5T EM- HEXANE AND BENZENE* ,/, OX , 26 ( "»H *),//» C X , 
1*KN OWN PO+NT— *» 2X,*P=*,F7.1»2X,*T0=*,F6.1»2X,*Xl=*,F7-4,2 v »*X2=*, r 
27.4,/) 

4 F0RMAT(1H1,/,50X,*SY2TEM- METHANOL AND BENZENE*,/ ,50X,2°( ’H*) ,/ t , c 
IX,* KNOWN P0+NT-*,2X,*P=*, F7.1,'>X,**0=*,F6. 1,2X,*X1=*, F7.4,2X, 

2 , F7 .4, / ) 

5 FOR MATtlHl, / ,50X,* SYSTEM- METHYLCYCLOPENTANE AND BENZEN C *, /, 5'?X,3° 

. 1 (1H*) ,//,5X,*KN0WN POI NT-* »2X,*P=*,F7.1,2X,*T 0=* ,F6.1,2 V , *Xl = *, f?. 

24, 2 X, *X2=*, F7.4, / ) 

6 FORMAT(lHl»/» 50X, * SYSTEM- ETHANOL AND METHYLCYCLOHEXANE*, ',50X,' 5 7 f 
11H* ),//, 5X,* KNOWN POINT- *,2X,*R=*, P 7.1»2X,*T0=*,F6.1,2X,* V 1.=*,F“'. A 
2,2X,*X2=*,F7.4,/) 

7 F0RMAT(F3.1 ,F4.2,F4.2, F4.2,F5.P) 

8 F0RMAT(F6.4,F6.4) 

9 FORMAT! F8.4,F8.4) 

11 FOR MAT(F5.3,F7.1,F6.2) 

12 • FOR MAT(F7.5,F3.3,F" T .3,F7.5,F8.' , ,F7.3) 

13 FOR MAT(F4.1 , F6 . 1 ) 

14 FOR MAT(5X,F5.1,3X,R7.1,5X, F7 .4 ,5X , F7 . 4 ) 



15 F0RFAT(5X,F9.4,5X,F9.4) 

16 FOR MAT (2X»F7 • ,6F 1 9 .4 ) 

17 FORMAT! 6F15. 4) 

19 FORMAT (F6.4,F6.4) 

R=8 0. 064760. 

READ7,F,RH0,RH0C,S,VQSTAR 
D06C0MP=1,6 
N=1 
I N=Q 

MN=G 
MA= C 

REAC8, XI »X2 
REAC9, VL1.VL2 
REAC12,A1,B1,C1,A2,B2,C2 
D0333 11=1,2 

REAC11,VSTAR(II),ESTAR(II) ,TC{ T I) 

333 CONTINUE 

I Ft MP. EQ. 1 ) MQ=6 
I FI FP » EO - 2 ) MQ=4 
I F{ MP.EQ.3)MQ=2 
I F{ MP . EQ .4 ) MQ = 2 
I F( NP. E0» 5 ) MQ=3 
I F( MP . EO. 6 ) MQ =2 
DQ2C0 M=1,MQ 
REA C13 ,T» P 

IF( MP • ECU 1 ) PR INTI , P » T, XI » X2 
I F{ MP.EG.2)PRINT2,P,T, XI, X2 
I F( MP.EQ.3)PRINT3,P,T,X1,X2 
I F{ PP.EG.4)PRINT4,P,T,X1,X2 
I F{ MP.EQ. 5)PRINT5,°,T, XI, X2 
I F ( PP.E0.6)PRINT6,P,T, XI, X2 
U ( M )=T+273-16 

Tl=Bl*(l./{ Al-ALOG { P )/ 2.303) )- r 1+2^3. 16 
T2=B2*(l./( A2-AL0G(P)/2.303) ) -r2+?"’3. 16 
D0400 1=1,2 

CALL V IR IAL ( U , RHO ,S,VOSTAR,VST*R,E c TAR,VC,TC,M,I) 

4 CO CONTINUE 

CALL VPC(T»A1,B1,C1 , A2 ,R2 , C2 , V*1 , V°2 ) 

GAM A1 = EXP { A LOG ( P/V D 1 ) + (F/(R*U( M ) ))*( (VC(1,1)-VL1 )*{P-VP’ ) M 
GAMA2=EXP(AL0G(P/VP2)+ (F/(R*U( M ) ) )*( ( VC( 2, 2 )-VL2 ) *( P-VP^ 5 1 ) 
PRI Ml 4, T, P , GAMA1 , GAM A 2 
IF! N— 1 ) 1 C , 1. 0 , 100 

10 CALL RK( GAMA1,GAMA2,X1,X2,MN,R^A,RKB,GE,R,AT,M) 

PRI Ml 5, RKA , RKR 

CALL VANLAR{IN,GAM*1,GAMA2,X1, V 2,VLA,VLB,GE,R,AT,M) 

PRI M15, VLA , VLB 

CALL SCAHAM ! GAMA 1 , GAMA2 , XI , , MA , SHA ,SHB , GE ,R, AT , M , VL’ , ''L2 3 

PRI Ml 5* SHA , SHB 
N=N + 1 
GOT0200 


GH 
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ICO RKA=RKA*U{M-1 )/U(M) 

RKB =RKB*U( M— 1 ) /U ( M ) 

VLA =VL A*U { M— 1 ) /U ( Ml 
VLB = VLB*U<M — 1 )/Ut Ml 
$HA=SHA*U(M-1 )/U(M) 

SHB = SH8*U ( M— 1 ) /U ( M) 

REAC19,X1,PQ 
222 X 2= 1 « - X 1 

RGMA1=EXP(X1*X?*(RKA+RKB*{ Xl-X* 5 ) ) +V 2* { RK A* { X2-X1 ) +RKB* { A . *X1* X<>1 . 

1 ) )) 

RGMA2 = EXR{X 1*X2*(RKA+RKB* ( XI -X*’ ) ) - *1* ( RK A* { X2— XI 5 +RKB* ( * . *X1 * a^-I . 

1 } )) 

C RGMA1 , RGMA2-ACT I VI T V COEFFICIENT 0 GIVEN BY R EDL I CH-K I STF R EHJATIG M 
C FDR COMPONENT ONE ANR TWO RESPECTIVELY 

VGM Al = EXP <VLA*X2*XV(VLA*Xl/VL°+X2 1**2) 

VGMA2=EXP(VLB*X1*XI /( X1 + VLB*XZ/VLA1**2) • 

C VGMA1 ,VGMA 2 -ACTIVITY COEFFICIENT 0 GIVEN BY VAN-LAAR EGUATTQ*' FOR 
C COMPONENT ONE AND TWO 

Z1=X1/(X1+X2*VL2/VI 1) 

Z2= ( X2*VL2/ VL! ) / ( XI +X 2 *VL2 / VLi I 
SGMA1=EXP(Z2*Z2*CSHA+2.*Z1*{ SH n *VLWL2-SHA) ) ) 
SGMA2=EXP(Z1*Z1*(SHB+2.*Z2*(SH**VL^/VL1-SHB) J ) 

C SGMA1 ,SGMA2-ACTIVITY COEFFICIENT 0 GIVEN RY SCAT CHARD-HAMMER F0UAT T 0 M 
C FOR COMPONENT ONE AND TWO 
VRG 1=VP1 *RGMAI 
VRG2=VP2*RGMA2 
VVG 1=VP1*VGMA l 
VVG2=VP2*VGMA2 
VSG 1= VP1*SGMAL 
VSG2=VP2*SGMA2 

PRI NT16» XI » RGMAlt RGMA2 » VGMA1 » VPMA2* SGMA1 » SGMA2 
PRIKTI7, VRGl,VRG2,VVGl,VVG2»VS r l,V°G2 
IF( Xl.GE.PQ) GCT020 r ' 

X1=X1+0.0005 
GOT G222 
2C0 CONTINUE 
600 CONTINUE 
STOP 
END 

C THIS SUBROUTINE CALCULATES THE V^POU 5 PRESSURES OF THE CO**P n NENTS 
SUBROUTINE VPC ( T , A1 , B?. , Cl , A2 , B*» C2, VP1 , VP2) 

VP1=EXP<2.303*< Al-fU/( T+Cl ) ) ) 

VP2=EXP<2.303*( A2-B2/I T+C2 )) ) 

RETURN 

END 

C THIS SUBROUTINE CALCULATES THE VTRIAl COEFFICIENTS 

SUBROUTINE VIRIAL ( U, RH 0, S , VOST* R » V °T AR , ESTAR, VC , TC , M , I ) 

DIMENSION U( 10) , VC 1 10, 10) , VSTA° { 5 ) ■» ESTAR ( 5) ,TC(5) 



J=2 

.77+0. 13 2 5* { TC (T) ) = ! £: > t (^.334) 

Rl= ( l.+S/(VST.*R(I ) ) )**'i 
R2= < 1 . + V0STAR/ (VST A -R { I ) ) ) **3 
R=I .987 

UNITS OF R - 3 .CAL » C. C . /GM.— MOVE 
RP1 = { RHODES TAR (I ) . 239) / ( R*(j { mi j 

RP2=S*RP1 11 

. . WSTAR=VSTAR(I )*«3 

I Ft ilLT.ZIGnT-?^ 1 "'^ 1 " 1 ' l * [EVp(RD ll-l-)-tR2-ltl|*(C X p(op^ ) _ 1 _ ) , 

999 r;tirn’ J , = < 3 ’” !VC( I_i ' 1-1 '.n i mvc ( i , n+vci i-i , i-D } 




END 


50 


THIS SUBROUTINE IS FHR rfdi rrw ncr- n 
SUBROUTINE RKIGAMAI.GW «'?' JS 

DIMENSION AT! 100) i,X * ' N tRKA » Rkb , GE, R, AT , M ) 

MN= FN+1 

IF! MN.GT.l JGOT050 
81= X2*X2 

Cl= (X1*X1*X2 + 5.*X1*X?*V?-X?) 

B2=X1*X1 " * 

C2= ( XI — XI *X2 * v .2~5 . *X1* XI $X 2 ) 

RKB = { !AL 0G{ GAMA 1)/B1)-{alqq/q^ »»*-,) /r?» w/n /n 
RKA = ( ALOGfGAMAl )/BT )-Cl*RKB/B V 2 )/lCl/B ^^B2) 

GOTC888 

aaa 

END 

C TH Bn R pnn??^ INE 15 FnR VAM-LAAR r QU A t I ON 

SUBROUTINE VANLAR ( Tfsj.GAMAl y i v 7 \n * u. 

DIMENSION AT('iOO) * * * £l, ^*-A»VL6»GF»R f AT f ! w ) 

in=in+i 

I F( IN.GT.l ) G0.T060 

.VLB = UL0G!GAMA2))*a*+(Ji4AL0rfrAMM 1! /! X1 * AL0G(GAMA1) 

G0TC777 'l.+(Xi*ALQGfGAM^l))/( X?*ALQG! GAMA2 )} ) **? 

f ^ 1X2* VLB/VLaI )*VLA^^* X2/ ^ EA ^ !iC ^ X ^ , ^ X ^ +X ^ :fC ^®‘^^^*-^^*^VLR i ! t X''/VLA)/(v] 
777 RETURN 
END 

C ' TH Rm»nn5?^I INE IS FnR S CATC HARD-HAMMER EQUATION 

DIMENSION" A“wo," ;4 " A1 ’ G4 " A2 ’'' l ’ X ’’ MA ’ SH4 ’ SHf, - GF - R -AT.''.>'Ll,VL-l 
Z 1= XI/ ( X1+X2* ( VL2/VL 1 ) ) 

Z2=.X2*( VL2/VL 1 ) / (X 7 +X2=M VL2/VL'* ) ) 
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MA=RA+1 

I F ( f'A. GT, 1 5 GOTO? 1 

SHB = ( A LOG ( GAM ' 1 ) / ( 7 2*Z2*( 1 „-2.*Zl ) )-ALOG { GAMA2 )/ f 2 . *Z I* 7 I*ZI*VL'V' r 
1L1) )/{ (2.*Zl*VlI/VL2/{l.-2.*Zi.1)-(fl.-2.*Z2)/(2.*Z2*VL2 'V* 7 ) } ) 

SHA =ALOG( GAMA 2 ) / ( 2.*Z1 *Z1*Z2*V1 2/VL1 )-SHB*{ ( 1.-2 .*Z2 5 / ( *> . -72* VL V'' 
1 LI ) ) 

GOT CD 5 

70 GE= R*AT( M )#( XIM Z2*Z2*( SHA +2 . * 7 1# ( VL1 /VL? )-SBA ) ) ) + v 2*{ Z1 *L' 
lSHB+2.*Z2*CSHA*{VL?/VLi)-SHB> ))) 

105 RET CRN 
END 

GENTRY 


7Z 
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♦ PROGRAM IH* 

* $ 


PREDICTION OF AZEOTROPIC BOILING TEMPERATURES AND COMPOSITI n N USING 
REDLICH-KlSTERtVAN-LAAR AND SC A T CHARD-HAMM ~R c QUATIO*’S 


C FOR CALCULATION OF VIRIAL COEFFI r IEN T S KREGLESKI S METHOD I c USED 

C LIQUID MOLAR VOLUMES ARE ASSUMED TO D E INDEPENDENT OF TFM D . AND PPE. 

C FOR VAPOR PRESSURE CALCULATION A^TOI^'E EQUATION IS USED 
C UNITS OF R (MM OF HG ) CC*C . ) / { r M-M n LE ) { K) 


DIMENSION U(1 ) ,VCH0,10) » VST A° ( 5 ) » EST AR ( 5 ) ,TC{5) ,AT{ j,Od ,UU{ 20 ) 
R=80. 06*760. 

2002 FOR MAT ( 100X , F 10.7 ) 

1 FOR MAT { / ,15X,*R-K CONSTANTS*, //,6X* *GA MA 1* , 7X , *G & MA 2* , 9 V , *PK A* , ^ D v 
1,*RKB*,10X,*T1*,11X,*T2*,/ ) 

2 FORMAT! /, 15 X,*VAN-LAAR CON STAN T $*, // , 6X , *GAMAi*, ''X, *GAM* 2*,9X , *VL * 
1*,1GX,*VL8* ,1 X,*T7*,11X,*T2*, /) 

5000 FORM AT ( 5X,3Fi ‘ .4, 5X , 3F13. 8 ) 

103 FORMAT! /,15X,*SCATCHARD-HAMMER CONSTANTS*, //, 6X, *G AM A 1*, 7 V , *GAM« 2* 
1,9X,*SHA*,10X,*SHB*,10X,*T1*,1’X,* T 2*, /> 

3 FORMAT! 1H1, // , 10X ,*C0MPAR I SION OF OBSERVED AND PREDICTE n " r EMPER AT* 1 
IRES BY R-K » VAN-LA AR AND SCATCHARD HAMMER EQUATIONS*,//) 

5 FORMAT!/, 15X,*F0R PTHA NOL* , / , 8 V , *V°1 *, 9X, *VC ! 1 , 1 ) * ,7X, *HE1 TA*,8Y,* 
1 VC! 1,2)*,/) 

a FOR MAT!/, 15 X, *FOR BENZENE* ,/,8 v ,*VPl*,9X,*VC!l»l )*,7X,* r 'EI TA*,8Y,* 
1 VC! 1,2)*,/) 

3 20 FOR MAT ( / , 15X , *FOR HEXANE* ,/, 8X,*VP1.*,9X, *VC ( 1 , 1 ) *,7X , *D e L T A* f 8X,*v 
ICQ ,2)*,/) 

1002 FOR MAT ( F6 .1 ) 

1003 FORMAT (100X,F7.4) 

321 FORMAT!/, 15X,*F0R METHANOL*, /, «X , *VP 1* , 9X, *VC { 1 , 1 ) *, 7X , *D r LTA* , RX , 

1*VC (1, 2)*,/ ) 

3 22 FORMAT! /, 15 X,*FOR METHYL CYCLQP C NT ANE *» /,8X,*VPl*,9X,*VCfl.l)*,7 Y ,* 
1DEL TA* ,8X ,*VC ! 1 ,2 ) *» / ) 

20 03 F0RMATC2X,F8.4,2X,F8.4,2X, F13.‘ 5 ,2X,F9-2,2X,F8.2,2X,F9.2, / I 

2004 FORMAT !/»2F14.4,3F76.4,F15.4) 

1005 F0RM>AT!F5.3 ,F5.3) 

6 FORMAT!/, 15X,*F0R C THY L ACETAT c *,/,8X,*VP2*,9X,*VC!2,2)*,~’X,*Dc)T' 
1*,8X,*VC!1, 2)*,/) 
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11 FORMAT!/, 15X,*F0R C YCLCHEX ANE *, / , 8 Y , *VP2* , 9X , *VC ( 2 , 2 ) *, "’X . *DF LT * * * 
18X, *VC (1,2)*,/) 

323 FOR MAT ( /,15X,*FOR ^ETH YLC YCLOH^XAN* - * , / ,8X , *VP 2*, 9X , *VC { **, ' , )*,7X,* r ' 
1ELTA*,8X,*VC{:,2)*,/) 

ICO, FOR MAT { / » 15X, *FOR 'BENZENE* , / , 3 v ,*V n 2*, 9X,*VC ( 2 *2 5 * ,7X , * n F’ T/i.*, o Y , * 
1VC( 1,2)*,/) 

7 FOR MAT(10X»F6.1,10 Y ,F6.1»10X»F<'.I, , OX,F6.],3X,F7.4,3X,F' 7 .', SX,? - *.'' 
1 ,3X ,F?.4»3X,F7.4,3 Y ,F7.4, / ) 

22 FOR MAT ( 1H1,/»S0X»*SYSTEM- ETHY1 AL r OHAL AND ETHYL AC ETA^E *, / , 5 " 
19 (1 H* ) »//,5X,*KN0W*' POINT-*,2X,*P=*,F7.1,2X,*TQ=*, F6.1, ‘ , X.*Xl = *, F” 

2.4, 2X» *X2=* » F7.4, / ) 

33 FORMAT! 1H1, /, 50X,*SYST EM- SEN Z r MF AND CYCLOHEXANE*, /, 50 Y , '’H 1H*) , ' 
1 / ,5 X,* KNOWN POINT—* , 2X ,*P = *, F7 . 1 , 2 v , *TO=*» F6. I ,2X , *X7= *. F". 4, 'X,* v 
22=*,F7.4,/) 

314 FOR MAT ( 1H1, / , SOX, * SYSTEM- HEXANE A M D BENZENE*, /, t '.jX,26r t H*),//» K X» 
1*KN CWN P0+NT-*,2X,*P=*,F7.1,2X,*TD=*,F6.1,2X,*Xl=*,F7.4.2 v ,*X2 = *, r ' 

27.4, /) 

316 FOR MAT ( 1H1, / , 50X,*RY2TEM- M!ETH A NOL AND BENZ £N '*, / , 50X , 2° ( ’ H* ) , / / , r 
IX,* KNOWN P0+NT-*,2X,*P=*,F7.1 ,'>X,* T 0=*,F6.1,2X,*Yi = *,F7.4, RX, *X*'=* 
2 , F7 .4 , / ) 

318 FOR MAT (1 HI, / , SOX, SYSTEM- METHYLCYCLOPENTANE AND BENZEN’"*, /, 50a, 3° 
1(1H*) » / / , 5X , *KNQWN POI NT-* , 2X , *P=*,F7.1,2X,*T 0=* ,F6.i,2 v ,*Xl=*, r 7. 
24,2X,*X2=*»F7.4,/) 

319 F0RMATUH1, /,50X,*EYSTEM- ETHANOL »ND METHYLC YCLOHEXANE*, ' , SOX ,’7 » 
11H* ),//,5X, *KNOWN POINT-*, 2X, * D =* , F7 . 1 , 2X , *TQ=*, F6.1»2Xt* v, =*,F' , . / 
2,2X,*X2=*,F7.4,/) 

77 F OR M A T ( / »5X »F‘ 7 .4»5 v »F7.4»5X,F9.4,5 Y ,F8.4,5X»F8.2,5X,F8. n » f ) 

99 FORMAT! F4.1,F*i. 3, F 5.3, F6.1 > • 

222 FOR MAT ( F3.1 , F4 .2 , F4. 2 , F4. 2 , F5 , F6. 2 , F6. 2 1 F6. 2 ) 

109 F0RMAT(F7.5,F8.3,F7.3, FT . 5 , F8 .*» , F7.3 ) 

666 F0RMAT(F5.3,F7.1,F6.2) 

66 FOR MAT { /, 5X , F 8 . 2 , 5X , F9 .2 , 5 X, F 8 . 2 , 5 Y ,F9 .2 , 5X , F9.4 , / ) 

4 FOR MAT(/,5X,*CHECK POINT*, 5X» *R=* , ^7. 1 , 2X , *T0=*, F6 . 1, 2X , * v l=* , F" 7 .* 
1 ,2X ,*X2=*,F7.4,/) 

104 FOR MAT ( /, 2X , * APPRO Y I MATE TEMPERATURE EVALUATED=*,F8.2,/1 
2000 FOR MAT!10X,F5.1»10X»F7.4,10X,F 7 .4) 

88 FOR MAT ( /» 9X, *OBS TFMIP* ,10X,*RK T E*, I2X, *VLTE*, 12X, *SHTE*, / ) 

AQ=C. 0000001 

VL1 =46.07/0.789 

VL2=88. 01/0. 901 

D06 00MP=1 ,6 

N=1 

MN=0 

IN=0 

MA=C 

MM= 1 

REAC222, F,RHO,RHOC,S,VOSTAR,RKT, VLT,SHT 

REAC109,A1,B1,C1,A?,B2,C2 

D03 33 11=1,2 

REAC666, VSTARUI) ,FSTAR ( 1 1) , TC ( II ) 
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333 CONTINUE 

I Ft y P . EG • 1 ) MG = 6 
I F{ yp .EG. 2 ) MG =4 
I F ( PP.EQ.3) MO=2 
I Ft MP.EQ.4) MQ = 3 
I Ft PP.EG.5)MQ = 2 
I F( PP.EC»6)MQ=2 
DC2 COM= 1 , MO 
DIF P IN = 10*) . 

JJ= ; 

NN= u 

REAC99,TC,X1, X2»P 
UU( P)=TO+273. 6 

Tl=Bl*(l,/{ Al-ALOGtP J/2.303) }-r 1+2^3.16 
T2=B2*(I./( A2-ALGG(P}/2.303) )- r 2+2^3. i 6 
U ( M ) = T 0 + 2 7 3 • 1 6 
I Ft N-I ) 10.10* '00 
10 DQ4C0 1=1,2 

CALL VIRIAL(U»RHO,S,VaSTAR,VSTAR,E*TAR,VC,TC,M,n 

4 CO CONTINUE ' ' ' 

DEL7A = 2.*VC ( 1 ,2 )-VC ( 1 , 1)- VC( 2 , *» ) 

T=TC 

CALL VPC{T,A1,B1»C1 , A2 ,B2 , C2 , V n l , V D 2 ) 

I Ft NP.EG.X) PRINT22»P»T0,X1 ,X2 
I Ft PP.EQ.l) G0T044 
l n l VL1 =78. /O. 879 
VL2=84. /0.779 

IF( MP.EG.2) PRINT?3,P,T0,Xi ,X2 
I Ft P P - EQ.2 ) GO 7044 
VLl=32./0.7924 
VL2=78. /O. 379 

I Ft yP.E0.3)PRINT316,P,T0,Xl,X2 
I Ft PP . EQ. 3 ) G0T044 
VL1 =84. /O . 750 
VL2 =78. /G. 879 

IFt yp.EG.4)PRINT3I. Q ,P,T0,Xl,X2 

I Ft PP.EG.4 JGOT044 
VL1 =46.07/0.789 
VL2 =98 ./0.769 

IFt yP.E0.5)PRINT3I9,P r T0,Xl,X2 
IFt yp . EG. 5 ) GOT044 
VLl=86./0. 65937 
VL2=78./0.879 

IFt yP.EQ.6) PRINT314, P, T0,X1,X2 
44 AA= ALOGt P/VP1 ) 

BB= F/ ( R*U( M ) ) 

CC= (VCt 1,1) -VL1 )*(P-VP1) 

DD=P*DELTA*X2*X2 
IFt yP.EO.l) PRINT5 
IFt yp.EQ.2) PRINTS 
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I F( f'P.E0.3)PRINT32 1 
I F ( VP. EQ.4) PRINT32? 

I F( MP.EQ.5) PRINTS 
I F ( PP.EQ.6) PRIMT32'' 

PRIM66tVPl.VC ( 1, i), DELTA, VC ( 1,2 ),VL1 
AAA=ALOG{ P/VP: ) 

BBB = F/I R*U( M) } 

CCC = I VC (2,?. )-VL2)*(P-VP2) 

DDD=P*DELTA*X'#X1 
I F( VP.EQ.DPRINT6 
I FI VP.EQ.2) PRINT11 
I Ft VP.EQ.6)PRIMT10 r '0 
I F C VP.EQ.3)PRINT!Ono 
I F( VP.FQ.4)PR1NT1GOO 
I F{ VP.EQ.5JPRI NT32' 5 

PRI M66,VP2,VCI2,2) , DELTA, VC (1,2) ,VL2 

GAMA1=EXP{AL0G{P/V°1)+ ( F/ ( R*U I M ) ) ) *( I VC 1 1 , 1 ) -VL3 } * ( P-V D 1 ' + P*DEi T " 
1*X2 *X2 ) ) 

GAM A2= EXP I A LOG I P/VP2J+ (F/ (R*U( M J ) )*( IVC(2,2) -VL ’ } * I P-V n 2 ' +P*DcLT * 
1*XI *Xl ) ) 

CAL L RK( GAMA1 ,GAMA2,X1 ,X2, MM, R*A, R*B,GE, R, AT » M ) 

PRI Ml 

PRI M77,GAMA1,GAMA’,RKA,RKB» T1,T2 

CALL VANLAR ( IN , GAMA 1 f GAMA2 , XI , v 2 , VL A , VLB , GE , R , AT , M ) 

PRI M2 

PRINT77,GAMA1,GAMA2,VLA, VLB,T;. ,T2 

CALL SC AHA M ( GAM AT , G AMA2 , XI , X** » MA , SHA , SHB , GE , R , AT , M , VL - * , ''L2 ) 

PRI NT IN 3 

PRI M77 , GAM A1,GAMA?»SHA>SHB»T1»T2 
N=N + 1 
PRI M3 
GOT C 200 

100 REA Cl 002, AT (M ) 

PRI M4» P,T0,X1,X2 
PRI M104, AT ( M ) 

PRI MSB 
ATI V— 1 ) =RKT 

555 RKA=RKA*ATIM-I)/ATCM) 

RKB=RKB*AT ( M-3 )/ATIM) 

ATI P-1 )=AT(M) 

RKT =AT ( M ) 

ZZ1 =RKA 
ZZ2 =RKB 

CALL RK{ GAMA 1 ,G AM A? ,X1 , X2 , MN , R* A , RKB ,G E, R , AT , M ) 

444 UIM )=ATIM) 

DOS 00 1=1,2 

CALL VIRIAL(U,RHO,S,VOSTAR,VST''R,E < 'TAR,VC,TC,M,I) 

SCO CONTINUE 

DEL TA=2.4VCI1,2}-VC(1» l)-VCt2,* , 3 
T=ATIM)-273 .16 
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CALL VPC{T,A1,B1,C1 , A2 ,B2 , C2 , V»1 , V°2 ) 

Yl= (VPl*GAMAr*Xl ) * C XP( fVLl-VC H » 1 ) ) * ( P-VPi ) / ( R*lM M ) ) )/P 
Y2= l.-Y], 

C Y1 » Y2— VA POUR PHASE COMPOSITION G c COMPONENT ONE ANO TWO 

VI* {(VLX-VC (1,1 ))*( P-VPI) )/(R*«T(MI) 

V2= (( VL2-VC (2,2) ) * f P-VP2 ) )/{R**T(M>) 

V=X 1*V1+X2*V2 

DEL SI =R*AT { M ) *ALOG ( P/VP1 ) / (Tl- A T( MI ) 

DELS2=R*AT( M) *ALOGfP/VP2) / (T2-«T(M) ) 

DEL S0=Xl*DELSI+X2*nELS2 

W*Xl*(ALOG( Yl/Xl) )+X2*(ALOG(Y2/X2) ) 

ATl=Tl*Xl*DELSl/DfcLSO 

AT2=T2*X2*DEL$2/DELSO 

AT3=GE/DELS0 

AT4=R*AT CM) *V/D£LSO 

AT5=R*AT(M)*W/DELSn 

ATI M ) =AT1+AT2-AT3-AT4+AT5 

BT=AT(M)-273.16 

DIF F=ABS { BT— TO ) 

DIF XY=A8S ( Y I-Xl ) 

AW= ABS ( W ) 

CHP = VP1#GAMA1*X1 + V D 2*GAMA2^X2 
EE=AT(M) 

FF= XI 
GG= X2 

OP* veil, 1) 

00= VC { 2,2) 

OR=VC( 1,2) 

0S= DELTA 
HH= Cl FXY 
PP= AW 

I F ( CIFF.LT.DIFMIN) GOTO 11 1 
I FC CIFXY.LE .0.0001 )G0T01004 
I F( AW.LE.AQ)G0T010n4 
Xl=Xl+0. CC05 
X 2=1. — XI 

I F{ Xl.GE.1.0 )G0TOlO04 
I FC JJ-1) 555,55,102 
10C4 I F( NN— 1 ) 300 ,700,1001 
111 DIF N IN=D IFF 

I FC JJ— 1)555,55,102 
3C0 RKTE=EE 

GMAR1=GAMA1 
GMA R2=GAMA2 
GER =GE 
VPR 1=VP1 
VPR 2= VP2 
CHP R=CHP 
REX 1 = F F 
REX 2=GG 
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)/AT(M) 
) /AT f M) 


0T= C P 
0U= CO 
OV=CR 

ox= cs 

OFX VR = HH 
A WR = PP 

REA Cl 002 » AT ( V } 

REA Cl 005 1 XI , X2 
AT( R-l ) = VLT 

difmn=ico. 

JJ= JJ + l 

MM= MM + 1 
NN= NN+1 

55 VLA =VLA*AT{ M-1 
VLB=VLB*AT( M- 
AT(P-1)=AT(M) 

VLT = AT( M ) 

Z Z1 = VLA 
ZZ2 = VL8 
CALL VANLARdf 
GOT C 444 
7C0 V LT E = E E 

GMA V1=GAMA1 
GHA V2=GAMA2 
GEV =GE 
VPV1=VP1 
VPV2=VP2 
CHP V=CHP 
VEX 1=FF 
VEX 2=GG 
PT=CP 
PU=CQ 
PV= CR 
PX=CS 
OFX YV=HH 
AWV = PP 

REA C1002 1 AT ( M ) 

REA 01005, XI , X? 

AT{ M-l )=SHT 

DIFMN=100. 

JJ= JJ+l 

nn=nn+i 

102 SHA=SHA4AT( M-1 )/ATf M) 

M-1 JXATCM) 
AT{ R-1)=AT(M> 

SHT = AT { M ) 

SHT E=£E 
GMA S1=GAMA1 
GMA S2 = GAMA2 
GES =GE 


GAMA1 , GAMA 2 »Xl, v 2»VLA t VLB»GE: 


t AT t M ) 
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VPS 1=VP1 
VPS 2= VP2 
CHP S=CHP 
SEX1=FF 
SEX 2 = GG 
QT=CP 
QU= CQ 
QV=CR 
QX= CS 
DFX YS=HH 
AWS = PP 
ZZI =SHA 
ZZ2 = SHB 

CAL L SCAHAM (GAMAT ,GAMA2 , XI » » MA , SH A ,SHB » GE » R» AT , M , VL’ » V L2 ) 

GOT C444 

10 01 PRI NT7,UU(M) ,RKTE,VLTE,SHTE,Rc v l,R^X2, VEX1 , VEX2 , SEX 1 , SE V 2 
C RKTE- AZEOTROPIC TEMPERATURE PRED T CTE n BY REDL ICH-K T STEP. E^U'TIOM 

C VLTE- AZEOTROPIC TEMPERATURE PREDICTS 0 BY VAN-LAAR EQUATION 

C SHTE- AZEOTROPIC TEMPERATURE PRED T CTE n BY SCATCHARD-HAMMER EQUATION 

C REX1, REX 2- AZEOTROPIC COMPOSITION PREDICTED BY REDL ICH-K I S T E° EON. 

C VEXl, VEX 2- AZEOTROPIC COMPOSITION PRE n I CTED BY VAN-LAAR EQ M . 

C SEX1, SEX2-AZE0TR0PIC COMPOSITION PREDICTED BY SCATCHARD-H * M M ER EQN. 

PRI NT5000,DFXYR,DFXYV,DFXYS,AWP, AHV, AWS 
C DFXYR-DI FFERENCE IN X AND Y OBT A T NED USING REDL ICH-K ISTER E P N. 

C OF XYV-D I FFERENCE IN X AND Y OBTAINED USING VAN-LAAR EQM. 

C DFXYS-D I FFERENCE IN X AND Y OBTAINED USING SCAT CHARD-HAMMER EQN. 

PRINT2004,RKA,RK8,OT,OU,OV,OX 
PRI NT2003»GMARl»GMAR2t GERt VP R1,VPR?*CH PR 
PRI NT2004,VLA,VLB,PT,PL»PV,PX 
PRI NT2003*GMAV1,GMAV2,GEV» VPV1 ,VPV?,CHPV 
PRI NT2004»SHA, SHB ,QT » QU » QV »QX 
PRI NT2003 ,GMAS1»GMAS2, GE S , VPS 1 , VPS2 , CHPS 
200 CONTINUE 
6C0 CONTINUE 
PRI NT9 

9 FO* MAT ( 1H1 ) 

STOP 

END 

C THIS SUBROUTINE CALCULATES THE V*POU* PRESSURES OF THE CO M P n NENTS 
SUBROUTINE VPC < T, A1 , B1 ,C1 , A2 ,B ? ,C 2, VP1 , VP2 ) 

VP1 =EXP{ 2.303* f Al-Bl/tT+Cl 3) > 

VP2=EXP( 2.303* (A2-B2/{T+C2m 

RETURN 

END 

C THIS SUBROUTINE CALCULATES THE V T R IAL COEFFICIENTS 

SUBROUTINE VI RIAL (U, RHO, S, VOST»R , VST AR , ESTAR , VC , TC , M, I ) 

DIMENSION UClOUVCf 10,10), VSTA”{ 5 ) , ESTAR (5) ,TC(5 } 

J=2 

B=D .77+0.1325* (TC< I ) )** (0.334) 

Rl= Q.+S/CVSTARU ) ) ) **3 



R2=(1.+VGSTAR/(VSTAR(I>) )**3 
R=: .987 

c UNITS OF R - G.CALf C.C./GM.-M 01 E,OK 
RP1 = ( RHO*ESTAR< I ) *0. 239) / ( R*U ( «) ) 

RP2 =S*RP 1 

WST AR=VSTAR ( I )**3 

VC( I, I >=R*WSTAR*( 1 .- ( R 1-1 . ) * { F V P (R n l )-l. )-(R2-Rl )*(£: XPC°P'’ )-l . ) J 
IF { I.LT.2)G0TC999 

vc( i-i,j )=( 2 .*vcu-i»i-i)*vcu,n )/(vc(i,i)+vcci-i,i-m 

999 RETLRN 

END 

C THIS SUBROUTINE IS F n R R EDL I CH— K T STEP EQUATION 

SUBROUTINE RK ( GAM AT , GAMA2 , XI » MN,RKA,RKB , GE , R, AT , M ) 

DIMENSION AT (' 00 ) 

MN= MN + 1 

IFC NN.GT.l ) G0T050 
B1=X2*X2 

Cl= (X1*X1*X2+5.*X1*X2*X2-X2) 

B2= X 1#X1 

C2= (X1-X1*X2*,X2— 5.*X1*XI*X2 ) 

RKB = ( ( ALOG( GAMA1) /f>l)-(ALOG(GA ,4 A2) /B2) ) / ( Cl/B 1-C?/ B2 ) 

RKA=< AL0G(GAMA1 >/Bl)-Cl*RKR/Bi 
GOT 08 88 

50 GE=R*AT< M)* (RKA*X1*X2 + RKB*X1*X' , *(X1-X2) ) 

GAMA1 = EXP(X1*X2*(RKA+RKB*( Xi-X*’) )+ v 2*( RKA#( X2-X1 )+RKB*(6.*Xl*X2-l 
1 ) )) 

GAMA2 = EXP(X1*X2*( RKA + RKB*( Xl-X* 5 ) ) -VI* ( RKA*( X2-X1 ) +RKB*( A. *XI*X2~1 
1 ) )) 

888 RETLRN 
END 

C THIS SUBROUTINE IS FOR VAN-LAAR r QU A T ION 

SUBROUTINE VANLAR ( I N , GAMA 1 , GAM*2 , X 1 , X2»VLA,VLB,G*=, R, AT, «) 
DIMENSION AT ( 100) 

I N= IN + I 

IFC IN. GT. 1 ) GOT060 

VLA= ( ALOG( GAMA1) ) *( 1.+ ( X2* ALOG f GAM* 2 )) / ( X1*ALGG( GAMAl ) ) 1**2 
VLB=< AL0G(GAMA2) )*(l.+(Xl*ALOGfGAMAI) ) / t X2*AL0G( GAMA2 )>) **2 
GOT 0777 

60 GE=R*AT(M)*(X'.+VLB*X2/VLA)*(Xl/(Xl+X2*VLB/VLA))*({VLB*X' , /' , LA)/( y l 

1X2*VLB/VLA) )*VLA 

GAMA1=EXP( VLA*X2*X2/ (VLA*X1/VL P +X2)**2 ) 

GAMA2=EXP( VLB*X1*X1 / (Xi+VL8*X2 'VLA) **2 ) 

777 RETLRN 

END 

C THIS SUBROUTINE IS FOR S CATC HARD— HAMPER EQUATION 

SUBROUTINE SC AHAM { GAMA 1 » GA MA2» V 1 tX’ , MA ,SHA, SHB ,GF » R» AT, »», 'M.1 , VL' ! ) 
DIMENSION AT ( 100 ) 

Z1=X1/(X1+X2*(VL2/VL1) ) 

Z2=X2*( VL2/VL:)/(XT+X2*(VL2/VL'' ) ) 

MA= NA + 1 

I F( NA.GT .1 ) G0 T 070 



SHR = ( AL0G(GAMA1 ) f { Z2*Z2*( 1 .~2,-*Zl 3 )“ALOG( GAMA2 ) / ( 2. *Z1* 7 1 * 12 * VL *V U 
1L1) )/( ( 2 . *Z 1 *VL 1/VL2/ { 1.-2:^1111 )-{ f 1 .-2. *Z2 ) / (2.*Z2*VL2 'Vi 1 ) } ) 

SHA = ALOG ( GA MA2 ) /( 2 . *Z 1*ZI# '22*V! 2 /Vt l )- SHB* { ( 1 .-2 . *Z2 ) / ( ** . *7 2* VL"/'' 
1 LI 1 ) 

G0TC105 

70 GE= R*AT(M)#(X1 : MZ2*Z2*(SHA+2.* 7 1 J M <:: HB 5 MVL1/VL2)-^HA) ) ) + v 2* ( Z1 *Z 7 * f 
1SH8+2.*Z2*( SHA*(VL' > /VL1)-SH8) ) 1 ) 

GAMA1=EXP( Z2*Z2*( SHA+2 .*Z1*( SH°*VL WL2-SHA ) 3 ) 

GAM A 2= EX P ( Z I *Z1 * ( SHB+2.*Z2*{ SH' , *VL**/VL1-SHB ) ) ) 

105 RETLRN 
END 


GENTRY 


-2 
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I'IS 


GOT Cl 05 

lSHB^Ii2i!I , sJj:i ( , ^;“*{SHA+2.*M.(5 H B*(v Ll /VL. 

gamai=exp( 7?* 77* f CN , L * )~shb ) n > 

o erf Z *"* Z 1 * ( SH8+? - '/VLP-SHA)) ) 


l ' ZI =!' Vj .*' 
VL? /V> •• ) } } 
*Z? ) 


/” 


)- c H4 ) } } + V 2 J ' ( l. *£' a * 


CFNTRY 


RETLRN 
END 


■*Z2*{Sf »*VL'*/VL1-SHB) > 


##a}rajc 


